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(JSTO)  of  the  Department  of  Defense  (DoD)  Chemical  and  Biological  Defense  (CBD)  Program. 
JSTO  is  also  the  Chemical/Biological  Technologies  (CB)  Directorate  in  the  Research  and 
Development  (RD)  Enterprise  of  the  Defense  Threat  Reduction  Agency  (DTRA).  Contract 
HDTRA1-10-C-0025  is  titled  Medical  Countermeasures  for  CBR  Agents. 

This  project  was  initiated  by  Ms.  Nancy  Nurthen,  of  the  Information  Systems  Capability 
Development  Division  (J9-CB)  and  continued  by  Dr.  Christopher  Kiley.  It  was  funded  under 
DTRA  Contract  Number  HDTRA1-10-C-0025  to  Gryphon  Scientific,  LLC,  with  subcontractor 
Applied  Research  Associates,  Inc.  (ARA).  The  target  application  for  the  product  of  this  contract 
is  under  the  auspices  of  the  Joint  Project  Manager  for  Information  Systems  (JPM  IS)  of  the  Joint 
Program  Executive  Office  for  Chemical  and  Biological  Defense  (JPEO-CBD). 
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Section  1. 


Executive  Summary 


A  model  for  the  treatment  of  plutonium-238/239  (Pu-238/239)  internal  contamination  by 
trisodium  calcium  and  trisodium  zinc  diethylenetriamine  pentaacetate  (DTP A)  was  developed  to 
evaluate  the  efficacy  of  various  treatment  regimens.  The  baseline  uptake  and  decorporation 
model  was  developed  by  adapting  existing  models  for  inhalation  exposure,  uptake/respiratory 
clearance,  and  physiologically-based  biokinetic  distribution  for  plutonium.  A  treatment  model 
was  developed  for  DTPA  decorporation  of  plutonium  that  simulates  the  chelation  of  DTPA  in 
the  body  and  its  enhanced  excretion.  A  radiation  dosimetry  model  for  calculating  doses  was  also 
developed  and  integrated  with  the  biokinetic  model  and  treatment  model.  Collectively,  the 
models  simulate  the  in  vivo  distribution  and  excretion  of  Pu  over  time  and  the  impact  of  DTPA 
treatment  (via  chelation  of  Pu  and  enhanced  excretion).  The  composite  Pu-238/239 
decorporation  model,  presented  in  this  work,  estimates  plutonium  deposition,  absorption, 
distribution,  retention,  excretion,  and  response  to  DTPA  treatment  in  adult  healthy  males  after  an 
inhalation  exposure.  The  model  further  estimates  the  acute  red  bone  marrow  and  lung  doses,  as 
well  as  the  whole  body  effective  dose  as  a  function  of  time.  The  long  term  carcinogenic  risk  from 
a  given  exposure  and  treatment  regime  can  be  estimated.  The  results  of  the  model  compare 
favorably  with  human  and  animal  data  and  provide  estimates  within  the  uncertainty  limits  of 
alternative  models  developed  by  the  International  Commission  on  Radiological  Protection 
(ICRP)  and  the  National  Council  on  Radiation  Protection  and  Measurements  (NCRP).  Outputs 
from  the  model  include  the  distribution  of  plutonium  deposited  in  the  respiratory  tract,  the 
amount  of  Pu  in  different  tissue  compartments  over  time  with  and  without  treatment,  Pu-238/239 
excretion  rates,  and  radiation  doses  to  critical  organs.  Calculations  from  the  model  may  be  used 
to  analyze  consequences  of  exposure  to  Pu-238/239,  in  particular  the  resulting  absorbed  radiation 
dose  from  the  exposure,  severity  of  pneumonitis,  and  risk  of  long-term  health  effects.  These  data 
can  be  compared  to  model  output  simulating  the  effect  of  treatment  based  on  different  initiation 
and  duration  times.  Comparison  of  radiation  dose  reduction  from  different  treatment  regimens 
can  be  used  to  analyze  the  efficacy  of  different  treatment  approaches.  The  model  may  facilitate 
interpretation  of  Pu-238/239  bioassay  data  and  aid  in  treatment  planning.  The  Pu-238/239 
decorporation  model  is  a  valuable  and  versatile  tool  for  assessing  the  effect  of  exposure  to  Pu- 
238/239  and  subsequent  DTPA  treatment. 
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Section  2. 


Introduction  and  Purpose 

Accurate  modeling  of  medical  countermeasure  efficacy  against  chemical,  biological,  and 
radiological  agents  (CBR  agents)  is  essential  to  understanding  the  vulnerability  of  our  war 
fighters  on  the  modem  battlefield.  In  helping  to  calculate  the  benefit  of  countermeasures, 
modeling  can  inform  data-driven  purchasing  decisions  and  logistical  tradeoffs.  In  this  study, 
Gryphon  Scientific  and  Applied  Research  Associates,  Inc.  (ARA)  have  developed  models  to 
predict  the  efficacy  of  medical  countermeasures  in  preventing  casualties  and  reducing  the 
severity  and  duration  of  illness  caused  by  a  variety  of  CB  agents  and  in  terms  of  radiation  dose 
reduction  and  reduced  long-term  health  effects  risk  for  R  agents. 

This  volume  (prepared  by  ARA)  is  one  of  ten,  describing  the  medical  countermeasure  models 
constructed  for  this  project.  This  volume  focuses  exclusively  on  the  modeling  approach, 
parameters,  and  calculations  used  for  the  medical  countermeasure  model  (MCM)  for  plutonium- 
238/239  (Pu-238/239)  internalized  contamination.  The  composite  model  includes  the  appropriate 
parameters  necessary  for  calculating  the  impact  of  trisodium  calcium  (Ca)  and  trisodium  zinc 
(Zn)  diethylenetriaminepentaacetate  (Ca-DTPA  and  Zn-DTPA,  respectively)  countermeasure 
treatment  of  Pu-238/239  internalized  contamination.  The  report  describes  how  the  calculations 
are  made  and  provides  examples  of  model  outputs  and  comparisons  to  data  and  other  model 
outputs  when  possible.  Other  volumes  describe  models  for  B.  anthracis  (volume  1), 
organophosphates  (volume  2),  cesium- 137  (volume  3),  F.  tularensis  (volume  4),  sulfur  mustard 
(volume  5),  americium-241  (volume  6),  Y.  pestis  (volume  7),  botulinum  toxin  (volume  8), 
plutonium-238/239  (volume  9)  and  vesicants  (volume  10,  an  expansion  on  volume  5). 

This  paper  presents  an  inhalation  exposure  model  for  calculating  the  deposition  fraction  of 
Pu-238/239  in  different  regions  of  the  respiratory  tract  from  a  given  air  concentration.  The  paper 
also  presents  an  uptake/clearance  model  for  Pu-238/239,  which  describes  the  rate  of  transfer  of 
Plutonium  (Pu)  from  the  respiratory  tract  to  the  systemic  circulation  or  the  gastrointestinal  (GI) 
tract.  Subsequently,  a  biokinetic  model  that  describes  the  systemic  distribution  of  Pu-238/239 
within  the  body,  a  decorporation  model  for  DTPA  binding  and  removal  of  Pu-238/239  in  the 
systemic  circulation,  and  a  radiation  dose  model  for  calculating  the  absorbed  radiation  dose  in 
the  whole  body  and  to  two  critical  organs,  the  red  bone  marrow  (RBM)  and  the  lungs,  is 
presented.  Each  model  approach  is  described  and  justified  along  with  the  assumptions  and  key 
parameters  that  are  implemented  in  each  model.  The  composite  model  connects  each  set  of 
calculations  and  is  collectively  used  to  calculate  the  acute  and  committed  radiation  doses,  which 
can  be  related  to  potential  adverse  health  effects.  The  radiation  dose,  with  and  without  DTPA 
treatment,  can  be  calculated  and  compared  to  determine  efficacy  of  treatment.  Different 
treatment  initiation  times  and  treatment  duration  times  can  be  evaluated  to  determine  the  impact 
of  different  treatment  scenarios  on  radiation  dose.  A  model  for  estimating  the  long  term 
carcinogenic  risk  from  a  given  intake  and  treatment  regime  is  included. 
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2.1  Summary  of  Deliverables 

Below  is  a  description  of  the  four  deliverables  assigned  for  each  agent  of  interest,  and  a 
description  of  what  is  included  in  each  deliverable.  The  third  numbered  item  (the  “MCM  Model” 
deliverable)  is  the  deliverable  presented  in  this  interim  report. 

1.  Modeling  Approach 

The  modeling  approach  deliverable  describes  each  of  the  parameters  that  we  anticipate 
including  in  our  MCM  model.  For  each  parameter,  a  description  of  the  approach  for 
developing  and  justifying  the  parameter  is  presented.  The  approach  is  based  on  prior 
knowledge  of  the  agent  and  on  the  general  types  of  data  available  for  each  agent,  but  specific 
citations  are  not  included,  as  it  is  a  preliminary  document. 

2.  Modeling  Parameters 

The  modeling  parameters  deliverable  defines  the  value  or  function  for  each  parameter  used  to 
develop  the  MCM  model.  Each  parameter  is  supported  with  a  description  of  the  rationale  for 
choosing  the  parameter,  including  any  scientific  evidence  used  in  parameter  development  or 
assumptions  that  were  made. 

3.  MCM  Model 

In  addition  to  the  information  already  developed  in  the  “Modeling  Parameters”  deliverable, 
the  “MCM  Model”  includes  a  description  of  the  model,  user  inputs,  the  model  calculations, 
and  the  model  outputs.  For  the  biological  agents,  the  report  is  an  accompaniment  to  a 
preliminary  implementation  of  the  MCM  model  built  in  Microsoft  Excel. 

4.  MCM  Model  in  a  Graphical  User  Interface  (GUI) 

The  MCM  models  will  be  implemented  in  their  final  version  in  a  graphical  user  interface. 
This  implementation  of  the  model  will  include  feedback  and  adjustments  made  after  review 
of  the  previous  deliverables. 
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Section  3. 


Background  and  Methodology 

3.1  Importance  of  Plutonium-238  and  -239 

Plutonium  (Pu)  is  primarily  a  man-made  actinide  element  (although  trace  amounts  of 
primordial  Pu  do  exist)  and  has  no  stable  isotopes.  Pu-239  is  a  fissile  material  that  is  used  to 
create  chain  reactions  for  use  in  nuclear  weapons  and  nuclear  reactors.  Release  of  Pu-239  into 
the  environment  has  occurred  from  nuclear  weapons  tests  and  accidents.  High  specific  activity 
Pu-238  can  generate  significant  amounts  of  heat  and  therefore  has  industrial  use  in  radioisotope 
thermoelectric  generators  (RTGs). 

Between  the  two  most  commonly  used  Pu  isotopes,  Pu-238  poses  the  greater  threat  for  use  in 
a  radiological  dispersal  device  (RDD)  because  it  has  been  used  in  industry  for  purposes  other 
than  nuclear  weapons  and  reactors  and  is  therefore  more  readily  accessible.  According  to  the 
National  Nuclear  Security  Administration  (NNSA),  the  NNSA  Offsite  Source  Recovery  Project 
had  recovered  over  2000  abandoned  or  unwanted  Pu-238  sources  worldwide,  totaling  over 
10,000  Curies  (Ci),  as  of  2007  (NRC  2008).  NNSA  efforts  indicate  that  a  significant  number  of 
previous  industrial  sources  exist;  many  of  which  have  been  abandoned  or  are  no  longer  wanted 
nationally,  and  more  importantly,  internationally.  NNSA  has  this  ongoing  activity  to  limit  such 
sources  existing  in  non-secure  environments,  where  they  could  potentially  be  obtained  and 
collected  for  malicious  purposes.  Nevertheless,  abandoned  and  unwanted  sources  not  yet 
collected  by  NNSA  remain  potential  RDD  threats. 

Pu-239  is  an  impractical  target  for  an  RDD  in  part  due  to  stringent  security  of  nuclear 
weapons  and  fissile  material.  The  Pu-239  associated  with  nuclear  reactors  is  in  a  matrix  with 
other  extremely  radioactive  material  such  that  it  renders  the  material  extremely  hazardous  and 
not  a  practical  source  for  use  in  an  RDD  that  could  be  encountered  “in  the  field”.  However,  Pu- 
239  could  be  an  inhalation  hazard  in  the  event  of  an  improvised  nuclear  device  that  does  not 
detonate  properly  but  results  in  the  spread  of  radioactive  material,  similar  to  the  weapon  accident 
that  occurred  in  Palomares,  Spain  in  1966  (Maydew  1997). 

The  primary  form  of  plutonium  used  in  industrial  applications  is  in  the  oxide  form, 
plutonium  dioxide  (PuCh),  which  is  considered  insoluble.  The  metal  oxide  is  used  in  RTGs 
(NRC  2008),  and  DOE  indicates  that  “Plutonium  dioxide  (PUO2)  may  now  be  the  most  important 
and  most  thoroughly  studied  of  all  plutonium  compounds.  Due  to  its  chemical  stability  and 
relative  inertness,  it  is  the  preferred  form  for  shipping  and  storing  plutonium  at  the  present  time” 
(DOE  2006).  Since  plutonium  readily  oxidizes  in  air  (DOE  2006)  and  in  an  RDD,  an  explosion 
would  likely  be  used  to  disperse  a  plutonium  monolith  (further  facilitating  oxidation);  any  form 
of  plutonium  used  would  likely  oxidize  into  PUO2.  Pu-238  decays  by  alpha  emission 
accompanied  by  a  low-energy  photon  to  uranium-234  with  a  half-life  of  87.7  years.  Pu-239 
decays  by  alpha  emission  followed  by  a  low-energy  photon  to  uranium-235  with  a  half-life  of 
over  24,110  years. 
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3.2  Plutonium-238/239  Decorporation 

The  pharmaceutical  countermeasures,  calcium  and  zinc  diethylenetriaminepentaacetate  (Ca- 
and  Zn-DTPA),  are  the  only  Food  and  Drug  Administration  (FDA)-approved  decorporation 
agents  for  Pu  and  the  actinides  americium  (Am)  and  curium  (Cm)  (FDA  2009).  Therefore, 

DTPA  (for  actinides)  and  Prussian  Blue  (for  Cs-137)  are  currently  part  of  the  Strategic  National 
Stockpile  (SNS)  for  which  significant  stores  of  the  agents  are  able  to  be  quickly  distributed  if  a 
significant  radiological  release  were  to  occur  (FDA  2009,  FDA  2011,  REMM  2012). 

It  is  worthwhile  to  note  that  the  US  Department  of  Health  and  Human  Services  Radiation 
Emergency  Medical  Management  website  for  “Managing  Internal  Contamination”  lists  the 
following  alternative  treatment  options:  “DFOA  [desferrioxamine],  EDTA  [ethylenediamine 
tetraacetic  acid],  and  DTPA  [diethylenetriaminepentaacetate]+DFOA”  (REMM  2011).  However, 
the  authors  explicitly  state  that  these  options  are  not  FDA  approved.  Their  source  reference 
(NCRP  Report  161,  NCRP  2008a),  used  extensively  in  the  development  of  the  approach 
presented  in  this  work,  indicates  that  DTPA  is  the  preferred  treatment  for  actinides,  such  as 
plutonium  and  americium.  The  main  body  of  that  guidance  only  provides  details  on  the  use  of 
DTPA  in  the  treatment  of  internalized  actinides.  However,  ethylenediamine  tetraacetic  acid 
(EDTA)  was  used  early  on  in  the  treatment  of  internalized  plutonium,  but  was  determined  to  be 
significantly  inferior  to  DTPA  (Norwood  1960).  Desferrioxamine  (DFOA)  was  investigated  in 
animal  models  for  improved  plutonium  chelation  with  limited  success  (Volf  1986,  1993). 
However,  DFOA  has  been  discussed  as  an  alternative  treatment  option  in  the  event  of  a  massive 
demand  for  actinide  chelation,  corresponding  to  depletion  of  DTPA  supplies.  Such  an  emergency 
“off-label”  use  of  DFOA  arises  from  the  fact  that  it  is  used  clinically  for  chelation  therapy  in  the 
event  of  iron  poison  (Chang  2011).  The  potential  exists  for  its  “off-label”  use  in  the  event  of  an 
emergency.  However,  very  limited  data  is  available  for  the  chelating  agents  EDTA  and  DFOA. 

The  structure  of  the  composite  model  presented  in  this  report  will  allow  for  implementation 
of  other  decorporation  models,  as  additional  data  becomes  available  for  new  decorporation 
agents.  This  utility  will  allow  for  the  comparison  of  different  decorporation  agents.  Because  the 
composite  model  uses  a  series  of  sub-models  and  is  built  on  biokinetic  data,  it  provides  the 
flexibility  to  more  easily  adjust  or  integrate  additional  parameters  applicable  for  novel 
decorporation  agents.  As  an  example,  there  is  an  effort  to  improve  current  chelation  methods 
through  formulation  of  an  orally  administered  form  of  DTPA  (Phan  2006a)  and  development  of 
novel  chelating  agents  is  ongoing  (Cassatt  2008).  As  data  becomes  available,  orally  administered 
DTPA  kinetics  can  easily  be  integrated  by  adding  parameters  for  absorption  of  DTPA  into  the 
circulation  where  it  can  chelate  Pu.  Likewise,  as  data  becomes  available  on  novel  decorporation 
agents,  terms  may  be  added  to  the  Pu  biokinetic  model  to  reflect  removal  from  relevant 
compartments.  The  biokinetic  model  foundation  provides  a  useful  platform  for  additional 
development.  The  other  components  of  the  model,  such  as  the  inhalation,  uptake,  and  dosimetry 
models,  also  remain  valid  and  provide  the  main  structure  for  additional  Pu  decorporation  models. 

DTPA,  shown  in  its  acid  form,  in  Figure  1,  is  a  chelating  agent  with  amino-  and  carboxylic 
acid  groups  that  can  coordinate  and  bind  actinides  with  relatively  high  affinity  and  enable 
enhanced  renal  excretion  of  these  radionuclides.  DTPA  can  be  administered  by  intravenous  drip 
or  direct  injection  into  a  vein  as  the  Ca  or  Zn  salt  in  a  5  milliliter  (mL)  hyperosmolar  solution. 
The  optimal  dosing  regimen  for  an  adult  is  1  gram  (g)  of  Ca-DTPA  on  the  initial  day  of 
treatment,  followed  by  1  g  Zn-DTPA  per  day  for  subsequent  treatments,  until  the  radionuclide 
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burden  is  adequately  reduced  (FDA  2004).  For  pregnant  or  lactating  women,  treatment  should 
begin  and  continue  only  with  Zn-DTPA  (NCRP  2008a). 


Figure  1.  Diethylene  Triamine  Pentaacetic  Acid,  DTPA 

Human  and  animal  data  indicate  that  DTPA  effectively  increases  the  rate  of  elimination  of 
plutonium,  which  then  in  turn  reduces  the  radiation  dose  absorbed  by  the  body.  While  DTPA 
reduces  the  absorbed  radiation  dose,  it  does  not  treat  the  effects  of  radiation  dose,  such  as 
hematopoietic  cell  damage  and  loss.  The  benefit  of  DTPA  treatment  depends  on  the  total  Pu- 
238/239  absorbed  into  the  body  and  the  time  to  DTPA  administration.  The  actinide  elements  can 
be  incorporated  into  bone  and  result  in  long  term  retention  of  the  nuclides.  Early  administration 
of  DTPA  can  minimize  the  radionuclide  deposition  in  the  bone  and  afford  greater  reduction  of 
radiation  dose  and  improved  countermeasure  efficacy. 

3.3  Composite  Model  Approach 

This  paper  presents  a  composite  model  for  determining  the  efficacy  of  DTPA  treatment  of 
internal  contamination  of  Pu-238  and  Pu-239.  Descriptions  of  the  parameters  used  in  each  of  the 
sub-model  components  are  provided.  The  sub-model  components  include: 

•  An  inhalation  exposure  model  for  calculating  the  fractions  of  plutonium  deposited  in  the 
three  primary  regions  of  the  respiratory  tract  from  a  given  air  exposure  concentration.  This 
model  provides  the  initial  body  burden  after  the  exposure. 

•  A  radionuclide-specific  uptake  model  for  absorption  of  plutonium  from  the  respiratory  tract 
to  the  circulatory  system  and  the  GI  tract. 

•  A  radionuclide-specific  biokinetic  model  for  calculating  the  distribution  of  plutonium  within 
the  body. 

•  A  radiation  dosimetry  model  for  calculating  the  absorbed  radiation  dose  to  the  critical  target 
organs  for  Pu-238  and  Pu-239.  This  model  is  used  for  estimating  acute  effects  and  the 
projected  whole  body  dose  over  50  years  for  long-term  effects. 

•  Dose  bands  assigning  projected  severity  levels  in  lung  doses  for  the  pulmonary  effects  of  a 
given  exposure. 

•  A  decorporation  model  of  DTPA  for  binding  and  removal  of  plutonium  in  the  systemic 
circulation. 

•  A  long-term  carcinogenic  risk  model  for  projecting  cancer  risk  from  a  given  intake  of  Pu-238 
or  Pu-239,  with  or  without  treatment. 
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Each  modeling  approach  is  discussed  and  justified  along  with  the  key  parameters  for  each 
model  set.  Collectively,  these  models  can  be  used  to  calculate  the  acute  radiation  dose  to  critical 
target  organs  for  determining  acute  effects,  as  well  as  the  long-term  committed  radiation  dose 
and  associated  carcinogenic  risk.  The  difference  in  the  radiation  dose  resulting  from  different 
DTPA  treatment  regimens,  including  varying  the  start  time  and  the  duration  of  administration, 
can  be  calculated  and  compared.  The  composite  modeling  approach  is  outlined  in  Figure  2.  The 
parameters  necessary  for  calculating  the  effect  of  DTPA  countermeasure  treatment  of 
internalized  Pu  are  further  outlined  for  each  model  component  in  the  next  section. 
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Figure  2.  Composite  Pu-238/239  Decorporation  Model  Components 
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Section  4. 


Plutonium-238/239  Decorporation  Model 


4.1  General  Assumptions 

4.1.1.  Exposure 

The  current  modeling  effort  addresses  the  battlefield  setting  in  which  an  RDD  is  deemed  the 
most  likely  threat  scenario.  Persons  involved  in  an  RDD  event  may  be  exposed  to  a  significant 
amount  of  radioactive  particulates,  which  would  potentially  result  in  skin  contamination  and 
inhalation  of  radioactive  nuclides.  Since  the  focus  of  the  current  work  is  the  development  of 
models  for  assessing  medical  countermeasures  for  internalized  radionuclides,  the  primary  route 
of  exposure  for  internalization  of  radionuclides  is  assumed  to  be  through  inhalation.  Some  of  the 
inhalation  exposure  will  result  in  the  eventual  transfer  of  material  to  the  gastrointestinal  tract  and 
is  addressed  through  the  inhalation  model.  Ingestion  is  not  currently  considered  a  stand-alone 
route  of  exposure  for  the  battlefield  scenario. 

The  amount  of  Pu  internalized  as  a  result  of  dermal  absorption  from  skin  contamination  is 
assumed  to  be  negligible.  Only  in  cases  of  heavy  contamination  of  a  soluble  radionuclide  would 
absorption  through  the  skin  result  in  any  significant  internalized  contamination  (NCRP  2008a). 
Plutonium  dioxide  (PUO2),  the  chemical  form  used  in  this  modeling  effort,  is  fairly  insoluble. 
Personnel  clothing  provides  significant  protection  from  skin  contamination,  and  external 
decontamination  can  be  readily  performed. 

External  exposure  to  radionuclides  from  an  RDD-type  scenario  can  potentially  result  in  a 
significant  contribution  to  radiation  dose  but  it  depends  on  the  specific  radionuclides  involved. 
For  the  radionuclides  Pu-238  and  Pu-239,  the  primary  contribution  to  radiation  dose  is  from  the 
alpha  emission  with  a  minor  contribution  from  low  energy  photons  or  electrons.  The  alpha 
particle  cannot  penetrate  past  the  epidermis  and  therefore  does  not  pose  a  significant  external 
hazard.  Furthermore,  any  countermeasure  treatment  for  the  external  dose  contribution  would  be 
targeted  toward  acute  radiation  syndrome  rather  than  decorporation.  As  stated  previously,  this 
particular  modeling  effort  addresses  internal  contamination  and  decorporation  treatment. 
Although  modeling  the  impact  of  treatment  on  the  dose  contributed  from  external  radionuclide 
sources  would  be  valuable  in  number  of  scenarios,  it  is  outside  of  the  scope  of  the  current  work. 

4.1.2.  Physical  and  Chemical  Form 

As  will  be  demonstrated  in  the  results  section,  the  physical  and  chemical  form  of  plutonium 
particles  are  critical  to  accurately  estimating  the  plutonium  distribution  and  resultant  radiation 
dose  from  an  inhalation  exposure.  Although  National  Planning  Scenario  estimates  particle  sizes 
of  1-150  microns  (pm)  result  from  an  RDD  (DHS  2004),  only  particles  less  than  10  pm  are 
considered  respirable  (ICRP  1994a).  Larger  particles,  up  to  30  pm,  may  be  inhaled  in  still  air  and 
even  larger  particles  may  be  inhaled  during  high  wind  speeds.  However,  very  large  particles 
settle  from  the  air  quickly  or  if  inhaled  are  lodged  in  the  nasal  region  (ICRP  1994a).  There  is 
insufficient  data  for  development  of  models  for  particle  sizes  greater  than  100  pm  (ACGIH  1999, 
ICRP  1994a).  The  inhalation  model  adopted  for  this  work  (as  described  in  section  4.2)  can 
potentially  model  mono-dispersed  particle  sizes  ranging  from  1  -  100  pm.  In  the  stand-alone 
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GUI,  only  1  and  5  micron  particle  sizes  were  implemented  to  minimize  the  computational 
demand.  These  particle  sizes  are  currently  implemented  in  other  radiological  models  (ICRP 
1994a  and  NCRP  2008b)  that  can  serve  as  reference  for  our  outputs.  ICRP  (1994b)  regards  1 
micron-sized  particles  to  represent  indoor  environmental  exposures  and  5  micron-sized  particles 
more  representative  of  occupation  exposures.  Recent  research  on  modeling  parameters  for  an 
RDD  that  included  particle  size  distributions  from  blast  was  reviewed.  Findings  indicated  that 
particle  sizes  from  1-3  pm  were  the  most  important  with  regard  to  eventual  inhalation  dose,  with 
a  peak  dose  at  2.6  pm  (Srimok  2006).  Therefore,  the  default  values  of  1  and  5  microns  for 
particles  implemented  in  our  work  should  be  sufficient  for  this  modeling  effort. 

Particles  will  deposit  in  different  areas  of  the  respiratory  tract;  larger  particles  will 
preferentially  deposit  in  the  upper  respiratory  tract  and  extrathoracic  region,  while  small  particles 
can  reach  the  lower  respiratory  tract  and  pulmonary  regions.  These  deposition  parameters  are 
estimated  using  the  respiratory  model  described  below.  This  approach  allows  the  user  to  examine 
the  influence  of  particle  size  on  deposition,  Pu  distribution,  and  dose. 

The  chemical  form  chosen  for  this  work,  Pu02  is  justified  by  the  text  presented  in  Section 
3.1:  “The  primary  form  of  plutonium  used  in  industrial  applications  is  in  the  oxide  form, 
plutonium  dioxide  (PUO2),  which  is  considered  insoluble.  The  metal  oxide  is  used  in  RTGs 
(NRC  2008),  and  DOE  indicates  that  “Plutonium  dioxide  (PUO2)  may  now  be  the  most  important 
and  most  thoroughly  studied  of  all  plutonium  compounds.  Due  to  its  chemical  stability  and 
relative  inertness,  it  is  the  preferred  form  for  shipping  and  storing  plutonium  at  the  present  time” 
(DOE  2006).  Since  plutonium  readily  oxidizes  in  air  (DOE  2006)  and  in  an  RDD,  an  explosion 
would  likely  be  used  to  disperse  a  plutonium  monolith  (further  facilitating  oxidation);  any  form 
of  plutonium  used  would  likely  oxidize  into  PUO2.” 

It  should  be  noted  that  while  PUO2  is  relatively  insoluble,  differences  in  solubility  of  Pu-238 
and  Pu-239  have  been  observed  and  reflected  in  the  ultimate  dissolution  and  uptake  of  Pu  from 
the  respiratory  tract  (NCRP  2008b).  Such  variances  depend  greatly  on  the  specific  activity  in  the 
particulate,  the  chemical  composition  of  the  particulate  as  a  whole  (i.e.  presence  of  carrier 
material)  and  size  and  shape  of  particulates.  Therefore,  the  outputs  generated  by  this  model 
should  be  considered  representative  of  Pu02  in  general. 

4.1.3.  Inhalation  Exposure 

As  described  above,  inhalation  exposure  was  considered  the  primary  route  for  plutonium 
intake  for  the  modeling  scenario  chosen  for  this  work.  In  the  case  of  PUO2,  which  has  quite 
limited  solubility,  the  precise  deposition  in  the  respiratory  tract  will  greatly  impact  the  rate  and 
amount  of  plutonium  reaching  the  systemic  circulation.  The  inhalation  model  described  below 
can  estimate  the  regional  deposition  location  of  plutonium  in  the  respiratory  tract  if  key 
variables,  such  as  particle  size  and  wind  speed  are  known.  Reasonable  default  parameters  are 
provided  when  specific  scenario  information  is  not  available.  We  compare  our  results  against 
other  established  models,  dose  coefficients,  or  case  studies,  by  choosing  parameters  that  most 
closely  resemble  the  specific  comparison  presented. 

4.1.4.  Ingestion  Exposure 

Given  the  RDD  scenario,  ingestion  is  not  considered  a  stand-alone  route  of  exposure.  A 
portion  of  the  plutonium  deposited  in  the  respiratory  tract  will  be  eventually  mechanically 
cleared  to  the  GI  tract  from  the  lungs  via  ciliary  action  or  swallowing.  PUO2  is  not  readily 
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absorbed  by  the  GI  tract  (0.001%,  ICRP  1993),  but  the  fraction  of  plutonium  reaching  the 
systemic  circulation  from  the  GI  tract  is  included. 

4.2  Inhalation  Exposure  Model:  Calculation  of  Plutonium  Body  Burden 

Medical  planners  determine  the  plutonium  body-burden  after  an  exposure  to  a  given  air 
concentration  of  plutonium  by  either  assumption  or  the  use  of  a  transport  and  dispersion  model. 
The  plutonium  body-burden  is  calculated  by  considering  the  inhaled  fraction  of  that  exposure 
and  the  amount  of  the  inhaled  fraction  that  gets  deposited  into  the  different  respiratory  regions. 
Both  of  these  parameters  depend  on  particle  size  and  can  provide  an  accurate  representation  of 
the  amount  of  plutonium  retained  from  the  exposure  (ICRP  1994a).  An  overview  of  the 
inhalation  model  used  to  calculate  the  plutonium  body  burden  is  provided  in  Figure  3.  The 
deposition  refers  to  respiratory  deposition  in  the  extrathoracic,  tracheobronchial,  and  pulmonary 
regions. 


Exposure  Inhaled  Fraction 


•Air 

•  Particle  size 

concentration 

•Wind  speed 

•Exposure 

duration 

•  Breathing  rate 

Intake  Deposition 

Fraction:  ET,  TB,  P 


m 

•  Inhalability 

m 

•Intake 

fraction 

•Particle  size 
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•Breathing  rate 

\ 

\ 

Figure  3.  Inhalation  Exposure  Model  Overview 


4.2.1.  Exposure 

The  exposure  (Exp)  is  a  necessary  input  for  the  inhalation  model  and  is  calculated  simply  by 
multiplying  the  air  concentration  of  plutonium  in  the  environment  (AirConc),  the  exposure 
duration  (ExpDur),  and  the  breathing  rate. 

Exp  =  AirConc  x  ExpDur  x  BreathingRate 

The  air  concentration  is  a  user-defined  input  parameter  in  units  of  activity.  In  the  final 
software  implementation,  variable  activity  inputs  will  be  possible  and  subsequently  converted,  as 
needed.  However,  for  the  basic  mathematical  construct  in  this  phase,  microcuries  (pCi)  per  unit 
volume  (liters)  are  used.  Also,  exposure  duration  and  breathing  rate  are  user-defined  inputs  in 
units  of  time  (minutes)  and  liters  per  minute  (L/min),  respectively.  The  user  may  want  to 
examine  different  levels  of  activity  and  the  impact  on  overall  calculations.  However,  for  this 
phase,  default  values  of  15  L/min  for  a  breathing  rate  at  rest  and  25  L/min  for  light  activity 
(Niinimaa  1981)  are  implemented. 

4.2.2.  Inhaled  Fraction 

Inhalability  is  a  measure  of  the  fractional  concentration  of  material  that  is  inhaled  relative  to 
the  concentration  in  the  ambient  air.  Inhalability  refers  to  the  intake  process  and  is  defined  as  a 
function  of  particle  size.  A  recent  review  (Millage  2010)  examined  the  data  from  inhalability 
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measurements  and  current  mathematical  models  of  inhalability.  Based  on  this  review,  the 
proposed  mathematical  model  provides  the  best  description  of  the  data  across  a  variety  of 
conditions,  while  minimizing  mathematical  inconsistencies.  The  inhalability  model,  which 
provides  the  inhaled  fraction,  is  a  function  of  both  particle  size  and  wind  speed;  wind  speed  can 
play  a  significant  role  in  how  much  particulate  material  can  be  inhaled. 

The  inhaled  fraction,  1(d),  is  determined  by  the  particle  size  id)  and  particle  concentration 

(C). 


1(d)  CinSpire(i(d)/Cami}ienf:(d) 


Inhalability  is  a  function  of  particle  size  and  ambient  wind  speed  and  is  best  described  by 
three  different  equations,  each  describing  different  wind  speed  regimes.  The  first  equation  is  for 
still  air,  i.e.  wind  speed  is  zero  (Millage  2010). 


I  still  (d)  —  //v 


(1  +  6.809  x  103  x  d~2-736) 


)  +  fo 


1.44 


(1  +  0.44  exp(0.0195d)) 


where  Ism  is  the  inhaled  fraction  in  still  air  and /n  is  the  nasal  fraction  as  defined  below. 


The  next  equation  is  for  winds  speeds  ( U)  greater  than  zero,  but  less  than  or  equal  to  4  meters 
per  second  (m/s)  (Millage  2010). 


if 

i 

1.44  l] 

.1  -  (1  +  exp(8.826  -  2.731n(d)))_ 

'  JO 

[l  +  0.44  exp(0.0195d)J  j 

+  J 


1-0.5  1 


-fe 


6  x  10_4d2-8  + 


-H  +  lx  10~SU2J5  exp(0.055d) 


U  2-75 

where:  /  =  717? 

The  final  equation  is  for  wind  speeds  greater  than  4  m/s  and  less  than  or  equal  to  9  m/s 
(Millage  2010). 

I^<lf^(d)  =  1  -  0.5  (l  -  (76xlQ-4d2.8  +  1))  +  1  x  10 ~5U2J5  exp(0.055 d) 


Note  that  data  is  not  available  for  estimating  the  inhaled  fractions  during  wind  speeds  greater 
than  9  m/s  and  the  model  is  not  valid  for  particle  sizes  greater  than  100  pm  in  aerodynamic 
diameter.  Therefore,  the  inhalability  function  is  limited  to  wind  speeds  of  9  m/s  or  less  and  to 
particles  sizes  of  100  pm  or  less.  Often  particles  larger  than  a  few  microns  are  ignored  because 
the  general  thought  is  that  they  will  not  be  inhaled;  however,  larger  particles  can  be  inhaled  and 
while  they  may  not  reach  the  pulmonary  region,  they  can  deposit  in  the  extrathoracic  region  and 
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be  absorbed  into  the  bloodstream.  As  an  example,  based  on  the  equation  shown  above,  in  a 
moderate  wind  environment  of  4  m/s  approximately  60%  of  22  pm  particles  will  be  inhaled. 

4.2. 2.1  Nasal  Fraction 

According  to  Niinimaa  1981,  the  majority  of  people  breathe  through  their  nose  until  their 
breathing  rate  reaches  a  threshold  level.  The  nasal  fraction,  fa,  defined  by  Niinimaa  is  illustrated 
in  Figure  4.  The  nasal  fraction  is  1  until  the  breathing  rate  reaches  greater  than  35  L/min,  which 
is  the  transition  point  to  heavier  activity.  The  oral  fraction,  fa  ,  is  simply  1  —  fN.  The  estimates 
established  by  Niinimaa  are  considered  to  be  representative  of  the  general  population  and  are  the 
best  estimates  outside  of  individual  data.  Note  that  this  approach  is  valid  only  for  normal 
augmentors  (i.e.  the  general  population).  In  the  final  software  implementation  of  the  model, 
values  will  be  implemented  in  a  look-up  table  for  simplicity. 


10  20  30  40  50  60  70  60  90 


Total  Respiratory  Minute  Volume  (I/m) 

Figure  4.  Nasal  Fraction  as  a  Function  of  Breathing  Rate 
(According  to  Niinimaa  1981,  reproduced  from  Millage  2010). 

4.22.2  Wind  Speed,  U 

If  wind  speed  conditions  are  known,  the  model  will  allow  a  user  to  enter  the  specified  value; 
otherwise  a  default  ambient  air  speed  value  of  4  m/s  will  be  used.  We  have  successfully 
accounted  for  varying  wind  speeds  in  our  implementation  of  the  particle  deposition  model  in  the 
Hazard  Prediction  &  Assessment  Capability  (HP AC),  so  the  final  implementation  of  the  model 
will  be  able  to  account  for  more  complex  weather  profiles. 

4.2.3.  Intake 

The  intake  is  simply  exposure  (in  activity  and  as  calculated  in  4.2.1)  of  radionuclide 
particulates  multiplied  by  the  inhaled  fraction  as  calculated  according  to  the  inhalability  function: 
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Intake  =  Exp  x  1(d) 

Where  Exp  is  the  exposure  and  1(d)  is  the  particle  size  and  wind  speed-dependent  inhaled 
fraction  which  valid  for  normal  augmentors  (i.e.  the  general  population). 

The  resultant  intake  is  the  radioactive  material  concentration  that  actually  enters  the  nasal  or 
oral  cavities.  It  does  not  take  into  account  the  fraction  of  the  material  that  is  actually  deposited  in 
the  respiratory  tract  or  the  fraction  that  is  breathed  out  again.  The  final  body-burden  will  be 
determined  by  the  deposition  in  the  respiratory  tract  as  described  below. 

4.2.4.  Deposition  Fraction 

The  deposition  fraction  is  the  fraction  of  the  inhaled  radioactive  material  that  is  deposited  in 
specific  regions  of  the  respiratory  tract.  We  used  the  Multiple-Path  Particle  Dosimetry  (MPPD) 
model  to  estimate  the  deposition  fraction.  MPPD  is  a  widely  used,  fast-running,  GUI-driven, 
Java-based  set  of  algorithms  that  can  calculate  the  deposition  and  retention  of  both  mono- 
dispersed  and  poly-dispersed  particulates  and  aerosol  droplets  in  human  respiratory  tracts  based 
on  user-defined  inhalation  rates  (see  Asgharian  2006a,  Asgharian  2006b,  and  Asgharian  2006c 
for  complete  details).  The  human  respiratory  model  includes  both  single-path  symmetrical 
calculations  as  well  as  several  multi-path  variations  of  limited-asymmetric,  asymmetric,  and 
stochastic  models.  The  limited-asymmetric  model  uses  a  5-lobed  model  with  subsequent 
symmetric  airways.  An  age-dependent  set  of  lung  morphologies  is  available.  The  model  can 
calculate  deposition  in  three  regions,  extrathoracic  (ET),  tracheobronchial  (TB)  and  pulmonary 
(P),  or  by  specific  airway  generation,  shown  in  Figure  5.  Clearance  pathways  and  time  frames 
are  significantly  different  in  each  region  and  since  Pu  is  relatively  insoluble,  the  estimated 
fractional  deposition  in  these  regions  can  significantly  impact  the  projected  clearance,  uptake, 
and  resultant  radiation  dose.  Therefore,  precise  deposition  calculations  are  critical  in  making 
accurate  estimations  in  the  Pu  model. 

Since  MPPD  does  not  rely  entirely  on  a  simple,  symmetrical  model,  such  as  that  used  in  the 
International  Commission  on  Radiological  Protection  (ICRP)  Publication  66  (1994a),  the  model 
can  account  for  specific  lobal  deposition  in  the  asymmetric  model  as  well  as  generation-by- 
generation  deposition  throughout  the  respiratory  tract.  Results  of  the  deposition  and  retention 
calculations  using  MPPD  have  been  favorably  compared  with  experimental  data  from  both  rats 
and  humans  (Raabe  1976,  Heyder  1986,  Asgharian  2006b).  Results  were  published  that  indicate 
the  asymmetrical  model  used  in  the  MPPD  is  effective  at  modeling  clearance  (Asgharian  2001). 
In  addition,  MPPD  includes  a  set  of  stochastically  generated  lung  models  that  can  be  used  to 
provide  an  estimate  of  the  uncertainty  associated  with  respiratory  tract  deposition. 
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Figure  5.  TheTthree  Functional  Regions  of  the  Respiratory  Tract  (Asgharian  2006a) 


As  an  example,  we  noted  previously  that  in  an  area  with  an  ambient  wind  of  4  m/s, 
approximately  60%  of  particles  of  22  pm  would  be  inhaled.  Using  MPPD,  we  calculated  that  of 
the  particles  that  were  inhaled  and  assuming  a  breathing  rate  of  15  L/min,  approximately  94% 
would  deposit  in  the  ET  region,  5.4%  would  deposit  in  the  TB  region,  0.0002%  would  deposit  in 
the  pulmonary  region,  and  the  remaining  0.59%  exhaled  without  depositing. 

In  much  of  the  published  data,  5  pm  particles  are  used  as  a  standard  for  inhalation  exposure 
scenarios.  In  cases  with  5  pm  particles,  100%  of  the  particles  would  be  inhaled.  Based  on  MPPD 
calculations,  we  calculate  that  approximately  31%  would  deposit  in  the  ET  region,  39%  would 
deposit  in  the  TB  region,  and  30%  would  deposit  in  the  P  region,  assuming  a  breathing  rate  of  20 
L/min.1 

Please  refer  to  the  references  used  in  this  section  for  a  comprehensive  description  of  the 
model  and  algorithms  included  in  MPPD.  MPPD  has  been  adopted  by  the  Environmental 
Protection  Agency  (EPA),  other  US  government  agencies,  and  the  international  risk  assessment 
community  to  precisely  estimate  respiratory  deposition  (TERA  2012,  Gangwal  2011,  Wichers 
2006,  RIVM  2002).  As  noted  above,  the  model  in  MPPD  has  been  documented  extensively  in 
peer-reviewed  publications  and  outputs  of  MPPD  have  been  transitioned  to  DTRA  programs  of 
record,  such  as  HPAC. 


1  The  breathing  rate  (20  L/min,  average  breathing  rate  for  a  worker)  and  particle  size  (5  micron)  were  chosen  for 
scenario  calculations  so  that  comparisons  could  be  made  with  other  published  data  that  use  these  parameters  (ICRP 
1994,  NCRP  2008b). 
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Since  MPPD  is  stand-alone  software,  a  more  comprehensive  look-up  table  using  MPPD  will 
be  implemented  in  the  final  software  developed  by  this  project,  to  ensure  that  all  components  of 
the  model  are  linked  together  in  a  seamless  fashion.  Users  will  not  be  required  to  run  MPPD 
separately  to  evaluate  variations  in  breathing  rate  and  particle  size  in  the  final  implementation.  A 
look-up  table  provides  a  simpler  and  faster  running  tool  than  implementation  of  the  MPPD 
engine  itself.  However,  for  illustrative  purposes,  Appendix  I  provides  an  example  of  the  outputs 
generated  by  MPPD.  Note  that  the  fractional  respiratory  deposition  depends  upon  both  particle 
size  and  breathing  rate. 

4.2.5.  Respiratory  Clearance  and  Systemic  Uptake 

The  size-dependent  deposition  profile  of  plutonium  particles,  along  with  other  key  attributes, 
will  impact  the  observed  clearance  and  systemic  uptake  of  plutonium.  For  example,  larger 
particles  are  less  likely  to  reach  the  pulmonary  region  in  which  absorption  into  blood 
predominately  occurs.  Deposition  in  the  extrathoracic  and  tracheobronchial  regions  of  the 
respiratory  tract  will  primarily  result  in  clearance  to  the  GI  tract  or  back  into  the  environment 
(via  nose  blowing,  coughing,  etc.).  Larger  particles  will  take  longer  to  dissolve,  which  is 
necessary  for  systemic  uptake,  thus  will  demonstrate  very  slow  absorption  into  the  bloodstream. 
Other  factors  which  impact  the  rate  of  uptake  include  the  chemical  form  of  plutonium, 
composition  of  the  matrix  material  of  the  particles,  and  presence  of  other  radionuclides 
(Menetrier  2005).  These  factors,  together  with  individual  variation,  contribute  to  the  large 
variability  in  the  data  available  on  the  respiratory  clearance  and  uptake  of  plutonium. 

Since  Pu02  is  relatively  insoluble,  a  fairly  small  amount  of  the  plutonium  deposited  in  the 
respiratory  tract  will  enter  the  systemic  circulation;  studies  indicate  that  plutonium  is  absorbed 
into  the  circulation  in  two  phases.  A  small  fraction,  less  than  1%,  of  the  deposition  is  absorbed 
within  about  a  day.  The  remaining  plutonium  is  mechanically  cleared  from  the  respiratory  tract, 
sequestered,  or  transformed  and  absorbed  into  the  blood,  with  a  half-time  on  the  order  of  years. 
The  rate  of  transfer  of  plutonium  out  of  the  respiratory  tract  can  vary  between  238Pu02  and 
239Pu02  as  demonstrated  by  bioassay  data  from  accidental  worker  exposures  (Hickman  2005  and 
James  2003)  and  from  animal  studies  (Bair  1973  and  Mewhinney  1983).  The  relatively  high 
specific  activity  of  Pu-238,  which  makes  it  amenable  for  use  as  a  RTG  can  cause  radiolytic 
fragmentation  of  Pu  particles,  increasing  solubility  and  facilitating  transfer  (NCRP  2008b). 
However,  due  to  limited  quantitative  data,  separate  values  for  Pu-238  are  not  currently 
implemented  in  our  model,  but  it  is  important  to  note  this  source  of  uncertainty. 

Considerable  evidence  exists  that  indicates  the  transfer  of  Pu  particles  is  slower  than  the 
generic  ICRP  (1994a)  respiratory  uptake/transfer  model  calculates  for  “insoluble”  radionuclides 
(Khokhryakov  2005,  Birchall  2010,  and  Davesne  2010).  The  impact  on  the  risk  to  pulmonary 
tissue  doses  (and  estimated  health  effects)  is  significantly  different,  depending  on  which  model 
parameters  implemented.  The  ICRP  (1994a)  parameters  for  respiratory  clearance  have  been 
critically  evaluated  since  additional  knowledge  has  been  obtained  and  improvements  to  the 
model  were  made  (Bailey  2007).  The  guidance  suggests  using  radionuclide  specific  parameters 
for  the  respiratory  tract  model  whenever  possible.  Plutonium  specific  parameters  were  estimated 
in  a  couple  of  studies,  since  the  ICRP  (1994a)  model  did  not  adequately  predict  values  obtained 
in  long-term  evaluation  of  exposed  plutonium  workers  (Khokhryakov  2005,  Birchall  2010). 
Likewise,  the  data  on  systemic  uptake  of  a  wide  variety  of  chemical  forms  of  plutonium 
particulates  from  human  and  animal  studies  have  been  critically  evaluated  and  specific  uptake 
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parameters  proposed  (Davesne  2010).  However,  each  of  the  referenced  studies  recommends 
different  parameters,  depending  on  the  specific  data  evaluated  or  the  approach  used. 

The  generic  structure  of  ICRP  respiratory  tract  model  for  respiratory  clearance,  as  illustrated 
in  Figure  6,  was  implemented  for  this  work.  The  clearance  parameters  are  listed  in  Table  1  below 
and  the  values,  shown  in  Figure  6,  refer  to  the  compartmental  transfer  rates  (in  day"1). 
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Figure  6.  ICRP  Respiratory  Tract  Model:  Mechanical  Clearance  (Bailey  2007)2 

The  ICRP  publication  provides  a  detailed  description  of  the  respiratory  tract  clearance  and 
uptake  model.  The  clearance  parameters  are  defined  in  Table  1,  with  the  fraction  of  regional 
deposition  in  each  compartment  provided.  To  correlate  with  outputs  from  the  MPPD  model 
described  previously,  the  ET  region  is  the  same,  the  TB  region  is  defined  as  bronchi  and 
bronchioles,  and  the  pulmonary  region  is  the  alveolar  interstitial.  Greater  detailed  outputs  from 
MPPD  can  be  obtained  in  terms  of  generations,  that  split  the  TB  region  into  bronchi  (generation 
1-9)  and  bronchioles  (generation  10-16).  These  outputs  from  MPPD  will  be  included  in  the  look¬ 
up  table  in  the  final  software  implementation. 


2  Bailey  (2007)  uses  a  prime  to  denote  the  portion  of  the  ET2  compartment  that  transitions  quickly  to  the  GI  tract. 
0.05%  of  the  initial  ET2  deposition  is  retained  in  the  wall  of  the  region  and  is  described  by  the  ETseq  compartment. 
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Table  1.  Respiratory  Compartments  for  Clearance  Calculations 


Compartment 

Symbol 

Fraction  of  deposition  in  each 
compartment 

Term  for 
calculation 

Extrathoracic  -  mainly  nasal 
deposition 

ETi 

0.15 

ET  1  frac 

Extrathoracic  -  lower,  cleared  to 
GI 

et2 

0.85*0.9995 

ET 2  frac 

Extrathoracic  -  sequestered, 
“bound” 

ETseq 

0.85*0.0005 

ET  seq  frac 

Bronchi  -  faster  clearance 

BBi 

0.993 -fs 

BBi  frac 

Bronchi  -  slower  clearance 

bb2 

fs 

BB2  frac 

Bronchi  -  sequestered,  “bound” 

BBseq 

0.007 

BBSeq  frac 

Bronchioles  -  faster  clearance 

bbi 

0.993 -fs 

bbi  frac 

Bronchioles  -  slower  clearance 

bb2 

fs 

bb2  frac 

Bronchioles  -  sequestered, 
“bound” 

bbseq 

0.007 

bbseq  frac 

Alveolar  interstitial  -  faster 
clearance 

AE 

0.3 

All  frac 

Alveolar  interstitial  -  moderate 
clearance 

ai2 

0.6 

AI2  frac 

Alveolar  interstitial  -  slower 
clearance 

ai3 

0.1 

AI3  frac 

Lymph  nodes 

LN 

LN  frac 

The  term  /v  refers  to  the  slow-cleared  fraction  in  the  “mechanical”  clearance  (i.e.  clearance 
mechanisms  besides  absorption/uptake  into  the  blood)  of  the  respiratory  tract,  which  is  particle 
size-dependent,  according  to  the  following  functions: 

fs  =  0.5  for  d  <  2.5  JpJx  pm 

-0.63  (d  1^— 2.5)  , - 

fs  =  0.5e  n*  for  d  >  2.5 ^fpjx  pm 

where  d  is  the  particle  size  diameter,  p  is  the  particle  density  (3  g/cm3  is  used  as  a  default 
value),  and  x  is  the  particle  shape  factor  (1.5  is  used  as  a  default  value;  ICRP  1994a). 

The  implementation  for  the  clearance  values  in  our  model  follow  that  of  the  ICRP  model, 
except  for  the  fraction  allocated  in  the  ET  region  for  clearance  to  the  environment  (via  nose 
blowing  or  coughing)  versus  clearance  to  the  GI  tract  or  sequestration.  ICRP  had  originally  used 
an  arbitrary  value  of  50%  clearance  to  the  environment;  however,  more  recent  data  indicates  that 
the  fraction  cleared  to  the  environment  is  no  more  than  about  15%  (Bailey  2007).  Therefore,  our 
model  assumes  a  15%  fraction  cleared  to  the  environment. 
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The  ICRP  model  assumes  that  absorption  into  the  bloodstream  occurs  in  all  regions,  besides 
the  ETi  compartment,  at  a  rate  specified  according  to  the  solubility  of  the  radionuclide  (soluble, 
moderately  soluble,  or  insoluble).  However,  the  revised  guidance  recommends  using 
radionuclide-specific  parameters  for  absorption/uptake  rates,  when  available  (Bailey  2007). 
Table  2  details  the  absorption/uptake  parameters  obtained  for  a  generic  insoluble  radionuclide  as 
well  as  from  different  studies  specifically  on  PuC>2.  A  couple  of  studies  have  proposed 
radionuclide-specific  parameters  for  the  clearance  model,  depicted  in  Figure  6,  together  with 
radionuclide-specific  parameters  for  systemic  uptake  (Khokhryakov  2005,  Birchall  2010). 
However,  these  studies  are  based  on  a  selected  case  or  experiment.  A  recent  review  by  Davesne 
2010  compiled  uptake  parameters  on  a  wide  variety  of  chemical  forms  of  plutonium  which 
includes  data  on  Pu-238,  -239,  and  -240  as  well  as  mixtures. 

Table  2.  Respiratory  Tract  Absorption  Values  for  Generic  Insoluble  Materials  and  PuC>2 


ICRP  1994; 
Insoluble 

Bailey  20073; 
Insoluble 

Birchall  2010; 
PuC>2 

Davesne  2010; 
P11O2 

fsr 

0.001 

0.01 

0.0026 

0.064 

sr  (day'1) 

100 

3 

1 

3 

ss  (day'1) 

0.0001 

0.00008 

0.0001 

0.0019 

where  fsr  is  the  percent  of  the  radionuclide  rapidly  absorbed  at  rate  sr  (day1)  and  the 
remaining  amount  of  radionuclide,  referred  to  as  fss,  is  equal  to  (l-/sy),  and  is  absorbed  at  a 
slower  rate,  ss  (day1). 

To  determine  an  appropriate  set  of  parameters  for  the  model  to  be  implemented  in  our  work, 
the  outputs  from  the  model,  using  each  of  the  different  absorption  parameters  listed  in  Table  2, 
were  compared  in  Figure  7. 


3  This  reference  uses  geometric  diameter  rather  than  aerodynamic  diameter  for  particles  in  their  estimate.  The 
authors  use  these  values,  there  may  be  differences 
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Pulmonary  Retention  of  Pu02as  Predicted  by 
Different  Absorption  Paramters 


- ICRP  1994 

Bailey  2007 

- Birchall  2010 

Davesne  2010 


Figure  7.  Pulmonary  Retention  of  PuC>2  as  Predicted  by  Different  Absorption  Parameters 


A  number  of  specific  parameters  impact  the  observed  behavior  of  plutonium  in  the 
respiratory  tract  such  as: 

•  Particle  size 

•  Exposure  parameters  such  as  wind  speed 

•  Inhalation  parameters  such  as  breathing  rate 

•  Chemical  form  of  the  matrix  material 

•  Physical  form 

•  Presence  of  other  radionuclides 

•  Individual  variability 

Considering  that  all  of  the  very  specific  parameters  (i.e.  particle  size,  presence  of  other 
radionuclides,  breathing  rate)  are  often  not  known-  and  for  the  purpose  of  scenario  testing-  a 
general  approach  that  encompasses  some  of  the  variability  is  warranted.  The  model  parameters 
provided  by  Bailey,  2007,  offer  a  more  general  response  as  compared  to  parameters  obtained  in 
the  specific  studies.  We  adopted  these  values  for  incorporation  into  our  model. 

The  portion  of  plutonium  transferred  to  the  GI  tract  from  the  respiratory  tract  via  mechanical 
clearance  (see  Figure  6  and  Appendix  II)  is  primarily  excreted  directly  through  fecal  elimination, 
although  a  small  fraction  of  the  plutonium  can  be  absorbed  into  the  systemic  circulation  from  the 
GI  tract.  Uptake  from  the  GI  tract  is  described  in  more  detail  below,  but  in  short,  studies  indicate 
that  plutonium  oxides  have  very  limited  absorption  from  the  GI  tract,  with  uptakes  ranging  from 
10"7  to  10"5.  A  GI  uptake  value  of  lxlO"5  has  been  adopted  by  the  ICRP  (1994b)  as  a  conservative 
estimate  for  insoluble  plutonium  uptake.  This  value  will  be  implemented  in  our  model  for 
consistency. 
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The  implementation  of  the  clearance  and  uptake  parameters  along  with  a  complete 
mathematical  description  of  the  model  is  provided  in  Appendix  II. 

4.2.6.  Uptake  from  the  GI  Tract 

The  portion  of  plutonium  transferred  to  the  GI  tract  is  primarily  excreted  directly  through 
fecal  elimination,  although  a  small  fraction  of  the  plutonium  can  be  absorbed  into  the  systemic 
circulation  from  the  GI  tract.  In  case  studies  and  experimental  studies,  plutonium  has  exhibited 
very  limited  absorption  from  the  GI  tract.  According  to  human  data  (from  Sellafield  releases 
(Hunt  1986)  and  from  a  review  of  GI  absorption  of  radionuclides  (Harrison  2001))  from 
individuals  who  consumed  seafood  contaminated  with  plutonium  and  americium,  0.00001-0.01% 
of  plutonium  absorption  from  the  GI  tract  was  observed.  Based  on  the  data  reviewed  for  PuCh, 
ICRP  (1994b)  adopted  a  conservative  value  of  0.001%  absorption  for  Pu02,  due  to  its  low 
solubility.  Therefore,  we  use  the  conservative  GI  uptake  parameter  value,  recommended  by 
ICRP  of  0.001%  in  our  model.  However,  it  is  important  to  note  that  Harrison,  et  al.  (2001) 
recommended  a  median  value  of  0.05%  for  generic  Pu.  A  recent  review  evaluating  uncertainties 
associated  with  Pu  dose  estimates  highlights  the  discrepancy  in  the  GI  absorption  parameters  to 
be  a  major  source  of  uncertainty  in  current  modeling  (Puncher  2012). 

4.3  Plutonium  Biokinetic  Model 

The  distribution  of  plutonium  throughout  the  body’s  tissues  is  calculated  according  to  Pu 
biokinetic  model  based  on  the  ICRP  (1993),  but  updated  and  modified  by  Leggett  (2005),  as 
shown  in  Figure  8.  The  modifications  made  by  Leggett  (2005)  include  refinements  in  the  tissue 
compartments  (such  as  blood  and  liver)  to  more  accurately  reflect  data  obtained  from  autopsy 
analyses  and  long-term  follow-up  of  occupationally-exposed  Pu  workers.  The  biokinetic  model 
implements  the  transfer  coefficients  detailed  in  Table  3. 

The  physiologically-based  model  developed  by  Leggett  was  chosen  because  it  was  the  most 
detailed  model  available  for  Pu.  The  physiologically-based  model  allows  more  detailed  and 
accurate  dosimetry  calculations  to  be  made  for  the  whole  body  and  enables  the  calculation  of 
doses  to  specific  tissues.  Specific  tissue  doses  are  critical  in  determining  acute  effects  from 
radiation  for  which  dose  to  critical  organs  must  be  examined.  Detailed  tissue  transport  data, 
obtained  from  the  biokinetic  model,  enables  the  calculation  of  the  amount  of  Pu  available  for 
removal  by  DTPA  treatment. 
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Figure  8.  Plutonium  Biokinetic  Distribution  Model  (Leggett  2005) 


The  plutonium  biokinetic  model  obtains  input  from  the  inhalation  model  to  ascertain  the 
amount  of  Pu  entering  the  bloodstream  from  the  respiratory  tract,  either  by  direct  absorption  into 
the  blood  or  indirect  respiratory  clearance  to  the  GI  tract  and  subsequent  absorption  (which 
contributes  minimally).  Due  to  the  low  solubility  of  PuC>2,  systemic  uptake  is  relatively  small. 
However,  significant  doses  to  internal  organs  may  still  be  obtained  from  the  portion  of  the  Pu 
absorbed,  particularly  from  Pu  immobilized  in  the  bone  compartment.  Therefore,  the  biokinetic 
model  is  required  for  calculating  those  doses  and  the  risk  of  subsequent  adverse  health  effects. 

4.3.1.  Transfer  Coefficients  for  Plutonium 

The  plutonium  biokinetic  model  allows  incomplete,  tissue-dependent  extraction  of  material 
during  passage  through  the  circulation  and  return  of  material  from  tissues  to  blood.  The  model 
performs  this  function  by  accounting  for  the  blood  flow  rate  to  different  tissues,  estimated  tissue 
retention  times,  and  subsequent  transfer  coefficients. 


Table  3.  Transfer  Coefficients  (Leggett  2005) 


From 

To 

Transfer  coefficient 

(d1) 

Blood  1 

Liver  1 

0.462 

Blood  1 

Cortical  surface 

0.08778 

Blood  1 

Cortical  volume 

0.00462 

Blood  1 

Trabecular  surface 

0.12474 
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From 

To 

Transfer  coefficient 
(d-1) 

Blood  1 

Trabecular  volume 

0.01386 

Blood  1 

Urinary  bladder  contents 

0.0154 

Blood  1 

Renal  tubules 

0.0077 

Blood  1 

Other  kidney  tissue 

0.000385 

Blood  1 

Upper  large  intestine  contents 

0.01155 

Blood  1 

Gonads 

0.0002695 

Blood  1 

Ovaries 

0.0000847 

Blood  1 

Soft  tissue  2 

0.018511 

Blood  1 

Soft  tissue  3 

0.0231 

Soft  tissue  1 

Blood  1 

0.099 

Blood  2 

Urinary  bladder  contents 

3.5 

Blood  2 

Blood  1 

67.55 

Blood  2 

Soft  tissue  1 

28.95 

Renal  tubules 

Urinary  bladder  contents 

0.017329 

Other  kidney 

Blood  2 

0.0001266 

Soft  tissue  2 

Blood  2 

0.001386 

Soft  tissue  3 

Blood  2 

0.0001266 

Liver  1 

Small  intestines 

0.0009242 

Liver  1 

Liver  2 

0.045286 

Liver  2 

Blood  2 

0.00152 

Liver  2 

Liver  3 

0.00038 

Liver  3 

Blood  2 

0.0001266 

Gonads 

Blood  2 

0.00038 

Ovaries 

Blood  2 

0.00038 

Cortical  surface 

Cortical  marrow 

0.0000821 

Cortical  surface 

Cortical  volume 

0.0000205 

Cortical  volume 

Cortical  marrow 

0.0000821 

Trabecular  surface 

Trabecular  marrow 

0.000493 

Trabecular  surface 

Trabecular  volume 

0.000123 

Trabecular  volume 

Trabecular  marrow 

0.000493 

Cortical  marrow 

Blood  2 

0.0076 

Trabecular  marrow 

Blood  2 

0.0076 

Urinary  bladder  contents 

Urine 

7 

GI  tract 

Feces 

0.5 
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The  transfer  coefficients  for  the  flow  of  Pu  from  blood  to  different  tissue  compartments  and 
its  return  to  the  circulation  were  mathematically  derived  from  the  fractional  blood  flow  to  the 
tissues  and  knowledge  concerning  the  tissue-specific  retention  of  Pu  in  the  compartments.  The 
knowledge  is  based  on  accidental,  worker,  and  even  experimental  human  exposure  case  studies 
and  autopsies,  as  well  as  detailed  experimental  animal  studies,  which  were  reviewed  by  ICRP 
(1993),  Leggett  (2005),  Durbin  (2006),  and  NCRP  (2008b). 

4.4  DTPA  Treatment  Model 

Since  DTPA  is  not  readily  absorbed  from  the  GI  tract  into  the  blood  stream  and  little  Pu  is 
secreted  into  the  GI  tract,  DTPA  must  be  administered  intravenously  for  effective  decorporation. 
Inhaled  DTPA  may  be  absorbed  up  to  20%  from  the  lung  into  the  circulation  and  is  a  potential 
mode  of  administration.  Based  on  data  from  animal  studies  and  one  human  tracer  study  using 
14carbon-labeled  DTPA,  the  chelate  uniformly  distributes  throughout  the  body  via  circulation.  It 
remains  in  the  extracellular  space  and  is  not  metabolized  to  any  significant  degree  (Volf  1978, 
Stather  1983).  A  basic  biokinetic  model  for  DTPA  has  been  proposed,  as  shown  in  Figure  9 
(Stather  1983,  Breustedt  2009  and  2010).  The  transfer  coefficients  proposed  for  this  model  are 
listed  in  Table  4  below. 


Figure  9.  DTPA  Biokinetic  Model  (Breustedt  2010) 

Breustedt  and  colleagues  through  the  European  consortium,  EURODOS,  have  tried  to  use  the 
biokinetic  model  to  develop  a  mathematical  model  for  DTPA  on  the  molecular  level.  To  date, 
these  efforts  have  had  moderate  success,  since  simulation  could  not  accurately  reproduce  data  for 
several  validation  cases  (Breustedt  2009  and  2010).  Due  to  these  limitations,  we  adopted  a 
simplified  approach  to  modeling  the  removal  of  plutonium  by  DTPA.  Since  DTPA  does  not 
concentrate  in  any  specific  tissues  or  transfer  to  intracellular  space,  we  assume  that  DTPA 
removes  Pu  by  chelating  Pu  primarily  in  the  circulation  (i.e.  the  blood  compartment  of  the  Pu 
biokinetic  model). 
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Table  4.  Transfer  Coefficients  for  the  DTPA  Biokinetic  Model  (Breustedt  2009) 


From 

To 

Transfer  coefficient  (d'1) 

Blood 

Interstitium 

145 

Interstitium 

Blood 

64 

Interstitium 

Lymph 

0.123 

Lymph 

Blood 

0.405 

Blood 

Urinary  bladder 

45.7 

The  process  of  decorporation  is  limited  by  the  rate  of  mobilization  from  the  lungs  to  the 
circulation  and  the  rate  of  release  of  Pu  back  into  the  circulation  from  tissues  that  have 
sequestered  it,  such  as  the  liver  and  bone.  While  the  Pu  biokinetic  model  calculates  the  amount 
of  Pu  in  the  circulation  at  any  given  time,  the  DTPA  biokinetic  model  developed  by  Stather 
(1983),  indicates  that  nearly  all  of  the  DTPA  injected  into  the  bloodstream  (99%)  clears  via  renal 
excretion  in  about  24  hours.  The  recommended  dose  of  DTPA  results  in  a  vast  excess  of  DTPA 
in  the  circulation,  as  compared  to  the  amount  of  Pu  in  the  circulation.  Also,  the  binding  affinity 
of  DTPA  for  Pu  is  quite  strong.  Both  of  these  qualities  assume  that  all  of  the  Pu  in  the  blood 
stream  will  be  bound  by  DTPA  during  its  circulation  and  is  subsequently  removed  by  renal 
clearance. 

4.4.1.  Route  of  Administration 

The  primary  route  of  administration  for  DTPA  is  intravenous  administration  and  is 
considered  a  fixed  parameter  in  the  model  (FDA  2004).  This  route  of  administration  will  quickly 
reach  an  equilibrium  concentration  of  DTPA  in  the  circulation  (Stather  1983). 

4.4.2.  Formulation,  Dose,  and  Time  Course 

Ca-  and  Zn-DTPA  are  formulated  in  a  5  ml  hyperosmolar  solution,  containing  1  gram  of  the 
respective  salt,  adjusted  to  physiological  pH  with  sodium  hydroxide  provided  in  sterile,  non- 
pyrogenic  ampoules,  suitable  for  intravenous  administration  (FDA  2004,  Hameln  2004).  This  is 
a  standard  formulation,  which  is  designed  as  a  fixed  parameter  in  the  treatment  model. 

Although  nebulized  DTPA  may  be  administered  with  an  inhaler,  there  is  insufficient  data  to 
adequately  model  this  route  of  administration.  The  evidence  that  is  available  does  not  indicate 
any  improved  chelation  for  Pu02  (Menetrier  2005).  In  addition,  nebulized  chelation  therapy  may 
exacerbate  asthma  (Hameln  2004);  therefore,  intravenous  injection  is  the  preferred  method  of 
treatment  and  is  treated  as  a  fixed  parameter  in  the  model. 

The  optimal  dosing  regimen  for  an  adult  is  1  g  of  Ca-DTPA  on  the  initial  day  of  treatment 
followed  by  1  g  Zn-DTPA  per  day  for  subsequent  treatments,  until  the  radionuclide  burden  is 
adequately  reduced  (FDA  2004).  The  transition  from  Ca-DTPA  to  Zn-DTPA  is  recommended 
due  to  toxicity  issues  and  is  briefly  discussed  later  is  this  report.  The  model  assumes  that  the 
recommended  dosing  regimen  is  followed,  although  the  mathematical  implementation  of  the 
decorporation  model,  as  described  below,  does  not  discriminate  between  Ca-  and  Zn-DTPA. 
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The  duration  of  treatment,  in  days,  is  an  input  variable  for  the  treatment  model.  The  effect  of 
treatment  duration  can  be  examined  with  the  model,  with  longer  treatment  courses  having  a 
greater  impact  on  total  radiation  dose. 

4.4.3.  Time  from  Exposure  to  Treatment 

The  time  from  exposure  to  initiation  of  treatment  is  an  input  variable  parameter  in  the  model. 
Using  the  plutonium  biokinetic  model  enables  the  dose  calculation  to  be  made  as  a  function  of 
Pu  concentration  in  the  body  integrated  over  time.  When  DTPA  treatment  begins,  Pu 
concentration  is  reduced  from  the  systemic  circulation  (the  blood  compartment).  Re-equilibration 
is  achieved  over  time  according  to  the  tissue  transfer  coefficients  in  the  biokinetic  model.  Since 
radiation  dose  is  calculated  based  on  the  Pu  concentration  as  a  function  of  time,  it  can  be 
calculated  seamlessly  using  the  model,  regardless  of  DTPA  treatment  start  time  or  duration.  The 
impact  of  time  from  exposure  to  initiation  of  treatment  can  be  examined  with  the  model.  The 
greatest  impact  on  dose  reduction  is  observed  when  treatment  begins  soon  after  exposure.  This  is 
particularly  important  for  Pu,  as  it  is  absorbed  into  the  bloodstream,  where  it  can  be  sequestered 
by  the  liver  and  incorporated  into  bone,  which  is  particularly  difficult  to  mobilize. 

4.4.4.  Plutonium  Decorporation 

The  binding  affinity  of  DTPA  for  Pu  is  quite  strong,  as  demonstrated  by  its  stability  constant 
of  23.4,  at  physiological  conditions  (Volf  1978,  Anderegg  2005).  The  stability  constant  is  the 
equilibrium  constant  for  the  formation  of  a  metal  complex  (a  chelate  bound  to  a  metal)  in  a 
solution  and  reflects  the  affinity  a  chelate  has  for  a  specific  metal.  For  comparison,  the  stability 
constant  of  DTPA  with  calcium  is  about  2  and  for  zinc  about  4  (Anderegg  2005).  Thus,  DTPA 
will  more  readily  exchange  zinc  or  calcium  ions,  allowing  for  the  opportunity  to  bind  to  Pu. 
DTPA  binding  of  Pu  is  stronger  than  with  other  competing  metal  ions  and  DTPA  does  not 
readily  release  Pu  once  it  is  bound. 

As  described  above,  the  recommended  dose  of  DTPA  results  in  an  initial  vast  excess  of 
DTPA  in  the  circulation  as  compared  to  the  amount  of  Pu  in  the  circulation.  Given  the  strong 
affinity  of  DTPA  for  Pu  and  its  excess,  we  can  assume  that  all  of  the  Pu  in  the  blood  stream  will 
be  bound  by  DTPA  during  its  circulation  and  subsequently  removed  by  renal  clearance. 

For  each  DTPA  treatment,  the  total  amount  of  Pu  entering  the  circulation  (blood 
compartment)  from  all  contributing  tissues  is  integrated  over  the  established  “effective” 
residence  time,  to  determine  the  total  Pu  bound  or  chelated  by  DTPA  for  each  treatment.  The 
amount  of  Pu  chelated  by  DTPA  is  deducted  from  the  total  body  burden  of  Pu  and  added  to  the 
amount  of  Pu  excreted.  The  Pu  remaining  in  the  tissue  compartments  is  then  re-equilibrated 
according  to  the  plutonium  biokinetic  model.  The  effective  residence  time  refers  to  the  period  of 
time  that  DTPA  can  efficiently  remove  Pu  from  the  body.  The  effective  residence  time  was 
developed  as  described  below. 

4.4.5.  Optimizing  DTPA  Effective  Residence  Time 

Different  effective  residence  times  for  Pu  in  the  circulation  (i.e.  blood  compartment)  were 
examined  in  the  treatment  model.  The  model  (described  above)  is  constructed  such  that  all  of  the 
Pu  entering  the  blood  compartment  during  the  effective  residence  time  for  DTPA,  after 
administration,  is  removed.  The  rate  of  Pu  chelated  from  the  circulation  during  DTPA  treatment 
is  the  sum  of  Pu  in  the  blood  that  would  normally  move  to  tissue  compartments,  and  the  amount 
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of  Pu  being  released  from  each  organ  back  into  the  circulation,  over  the  specified  effective 
residence  time  for  each  treatment: 


df-  chelated)  /  iP ublood  *  TC blood^orgarii  P^organi  *  TC organs  blood) 

where  Pucheiated  is  the  amount  of  Pu  chelated  by  DTPA,  Pubiood  is  the  amount  of  Pu  in  blood, 
Puorgani  is  the  amount  of  Pu  in  each  organ  (or  tissue)  that  will  transfer  within  the  blood  (see 
Figure  8),  and  TC  is  the  transfer  coefficient  (in  units  per  day;  Table  3)  for  the  respective  flow  of 
Pu.  The  amount  of  Pu  chelated  by  DTPA  is  added  to  the  amount  of  Pu  excreted  in  the  urine. 


The  transfer  of  Pu  in  and  out  of  the  blood  during  untreated  periods  is  represented  as: 
d 


a  v-1 

—  (Pu^;00(j)  =  y  (—fid blood  *  PCl 


blood->organi  Puorgarii  *  PCorgani->blood ) 


The  transfer  of  Pu  in  the  blood  compartment  during  treatment  changes  to: 

d 
dt 


'  (PUfriood)  ^^(  PV-biood  *  PCl 


hlood-^organi 


) 


Plutonium  retention  under  DTPA  treatment  was  examined  assuming  a  maximum  residence 
time  of  24  hours,  based  on  the  biokinetic  model  for  DTPA  (Stather  1983).  Data  was  generated  on 
effective  residence  times  for  DTPA  ranging  from  2  to  24  hours,  as  shown  in  Figure  10.  For  this 
scenario,  treatment  was  initiated  one  day  after  exposure  to  5  pm  particles.  Treatment  was 
continuous  for  30  days.  The  simulation  shows  how  effective  continuous  treatment  can  be  when 
initiated  early.  In  this  scenario,  Pu  is  removed  nearly  as  soon  as  it  enters  the  bloodstream.  Since  a 
slower  absorption  of  Pu  into  the  blood  occurs  as  particle  size  increases,  DTPA  treatment  can  be 
very  effective  if  administered  early  (before  Pu  becomes  lodged  in  the  bone  and  liver)  and 
continuous.  Since  Pu  moves  into  the  bloodstream  relatively  slowly,  different  effective  residence 
times  for  Pu  do  not  have  a  large  impact  on  the  amount  of  Pu  that  can  be  removed  with  each 
treatment. 
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Effect  of  Different  Residence  Times  of  DTPA  on 


Figure  10.  Effect  of  Different  Effective  Residence  Times  of  DTPA  in  the  Blood 

Using  the  plutonium  biokinetic  and  decorporation  model,  the  maximum  effective  residence 
time  for  24  hours  provides  an  initial  increase  of  approximately  7 1  times  the  daily  urinary 
excretion  of  Pu  during  treatment.  This  effectiveness  of  the  chelation  treatment  using  DTPA  is 
evaluated  by  considering  the  excretion  enhancement  factor  (EEF).  DTPA  has  been  reported  to 
increase  urinary  excretion  of  Pu  with  an  EEF  of  5  to  150  with  DTPA  treatment  (Jech  1972,  Jolly 
1972,  Schofield  1974,  Hall  1978,  Volf,  1978,  Wood  2000,  Menetrier  2005,  Grappin  2007, 
Bertelli  2010).  The  effectiveness  of  DTPA  in  six  patients  with  inhalation  exposure  of  239Pu 
(chemical  form  not  reported)  was  examined  five  years  after  exposure  with  intermittent  treatment. 
The  EEFs  ranged  from  45-120,  with  decreasing  effectiveness  over  time.  The  output  of  the  DTPA 
decorporation  model  (EEF-71)  falls  well  within  these  reported  values  for  enhanced  urinary 
excretion. 

Studies  on  the  treatment  of  Pu  and  Am  with  DTPA  in  animal  models  have  indicated  that 
DTPA  can  effectively  remove  actinides  from  tissues  in  the  extracellular/interstitial  space  and 
possibly  competitively  bind  nuclides  from  cellular  surfaces  (Phan  2005,  2006a,  2006b;  Fritsch 
2008,  2009,  2010).  Based  on  these  studies’  findings,  researchers  hypothesized  that  a  small 
amount  of  DTPA  can  actually  enter  intracellular  spaces  and  bind  to  and  remove  actinides.  This 
intracellular  attachment  mechanism  is  plausible  and  can  explain  the  observation  of  a  continued 
elevated  excretion  for  a  period  of  time  after  treatment  ends.  Although  the  vast  majority  of  DTPA 
clears  the  blood  in  24  hours  (99.9%),  a  small  amount  of  DTPA  clears  more  slowly  (Stather  1983) 
and  therefore  this  hypothesized  mechanism  needs  to  be  validated.  Sufficient  evidence  exists  to 
support  interstitial  removal  of  actinides  via  DTPA  treatment;  however,  all  of  the  reported  studies 
were  conducted  on  soluble  species  of  Pu  and  Am. 

To  account  for  an  elevated  increase  in  Am  excretion  after  DTPA  treatment  via  chelation  in 
the  interstitial  space,  removal  of  interstitial  Am  was  approximated  by  using  existing  transfer 
coefficients  from  the  different  tissue  compartments  and  estimating  an  interstitial  removal  fraction 
(IRF).  However,  in  the  plutonium  model  provided  here,  model  outputs  provide  an  EEF  (71)  of 
the  same  magnitude  as  observed  in  human  case  studies,  indicated  above.  Serandour  et  al.  (2008) 
proposed  that  the  reason  for  the  differential  behavior  of  Am,  compared  to  Pu,  indicates  a  lower 
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affinity  of  DTPA  for  Am  as  compared  to  Pu.  However,  the  composite  data  suggests  that  a  larger 
intracellular  contribution  to  decorporation  of  Am  excretion  under  DTPA  treatment,  as  compared 
to  Pu. 

It  is  worth  noting  that  a  decline  in  the  enhanced  excretion  is  observed  over  time,  as  tissue 
reservoirs  with  faster  transfer  coefficients  (reflecting  tissues  that  bind  plutonium  less  tightly)  are 
depleted.  In  turn,  the  excretion  enhancement  obtained  from  DTPA  declines.  Case  studies  show 
that  upon  cessation  of  treatment,  whole-body  re-distribution  of  Pu  from  slower  tissue  transfer 
compartments  occurs  and  other  soft  tissue  levels  rise  upon  re-equilibration.  Therefore, 
intermittent  treatment  with  DTPA  is  warranted  for  overall  effectiveness.  The  current  model  code 
is  constrained  by  the  software  used,  (see  section  4.1)  such  that  scenarios  with  intermittent 
treatments  are  not  easily  evaluated.  Future  implementation  may  include  this  feature. 

4.4.6.  Toxicity/Side  Effects 

Treatment  with  DTPA  can  deplete  endogenous  metals,  which  is  why  the  Zn  salt  is  preferred 
for  long-term  treatment.  However,  other  endogenous  metals  may  be  depleted  if  treatment  spans 
several  months.  During  prolonged  treatments,  serum  levels  of  essential  metals  should  be 
monitored  including  the  use  of  mineral  replacements,  as  needed.  Otherwise,  DTPA  has  not  been 
reported  to  have  any  significant  adverse  effects.  Approximately  6%  of  patients  in  the  registry  of 
DTPA-use  reported  side  effects,  such  as  headache  and  lightheadedness.  The  use  of  nebulized 
DTPA  for  inhalation  dosing  reported  patients  experiencing  coughing  and  wheezing  and 
exacerbated  asthma  (FDA  2004,  Hameln  2004).  Patients  with  renal  impairment  do  not  need  dose 
adjustment;  however,  if  renal  insufficiency  or  renal  failure  is  present,  Ca-DTPA  should  not  be 
used.  If  patients  with  renal  impairment  are  heavily  contaminated,  using  dialysis  increases  the  rate 
of  elimination.  Dialysis  fluid  will  become  radioactive  and  appropriate  radiation  safety  measures 
should  be  instituted  (NCRP  2008a). 

Although  long-term  treatment  with  Zn-DTPA  in  healthy  adults  does  not  cause  any  serious 
side  effects,  the  route  of  administration  is  through  intravenous  injection.  This  is,  to  some  degree, 
invasive  under  continuous  treatment  and  must  be  considered  in  the  overall  regimen  adopted. 

4.5  Radiation  Dose  Model 

The  radiation  dose  delivered  to  the  critical  target  organs  (RBM  and  lungs)  for  acute  effects 
must  account  for  the  Pu  burden  in  each  of  the  tissues.  The  amount  of  energy  absorbed  in  the 
critical  target  organ  must  account  for  the  energy  emissions  from  all  sources.  The  method  used  for 
this  calculation  is  sometimes  referred  to  as  the  Medical  Internal  Radiation  Dose  (MIRD)  method. 
The  method  utilizes  tabularized  data  that  are  computational  estimates  of  the  amount  of  energy 
deposited  in  a  specific  organ,  based  on  the  amount  of  energy  emitted  from  radioactive  material 
within  another  organ  (Snyder  1978).  Due  to  the  complex  nature  of  the  human  anatomy,  these 
specific  absorbed  fractions  (SAFs)  in  each  of  the  tissues,  require  detailed  radiation  transport 
calculations,  as  described  below.  The  dose  to  specific  target  organs  requires  an  integrated 
calculation,  based  on  the  dose  rate  as  a  function  of  time,  resulting  from  the  time-dependent  Pu- 
238  or  Pu-239  body-burden.  Since  Pu-238  and  Pu-239  both  emit  low-energy  photons  and 
electrons  (with  a  given  percent  per  disintegration)  and  an  a-particle,  the  dose  to  the  red  bone 
marrow  and  lungs  will  include  contributions  from  photons  emitted  from  other  organs.  Also, 
plutonium  within  the  RBM  and  lungs  will  further  emit  photons,  electrons,  and  a-particles. 
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The  local  energy  deposition  from  the  a-particle  is  modified  by  a  Quality  Factor  (QF)  of  20  to 
convert  to  dose  equivalent.  The  high  linear  energy  transfer  (LET)  of  the  a-particle  causes  a 
higher  level  of  localized  damage  than  a  low  LET  energy  deposition  from  photons.  To  quantify 
the  increase  in  potential  biological  damage,  multiply  the  energy  deposition  by  a  QF.  For  a- 
particles  in  the  5  mega-electron  Volts  (MeV)  range,  the  accepted  QF  is  20  (ICRP  1991). 

In  addition  to  the  organ- specific  dose  calculation,  we  present  a  method  for  approximating  the 
50-year  committed  effective  dose  equivalent  (CEDE).  The  CEDE  is  the  dose  calculation  that  is 
commonly  used  for  radiation  protection  associated  with  internal  hazards  and  it  is  the  unit  upon 
which  most  long-term  complications  or  cancer  incidence  rates  are  based.  Acute  dose  calculations 
are  generally  considered  to  be  associated  with  the  dose  accumulated  over  the  first  thirty  days;  the 
majority  of  the  total  dose  accumulated  will  occur  during  this  time  period.  The  body-burden,  and 
the  associated  dose  accumulation  rate,  will  decrease  with  time  as  a  result  of  radioactive  decay, 
natural  elimination,  and  accelerated  elimination  resulting  from  DTPA  treatment. 

4.5.1.  Organ  Specific  Acute  Doses 

The  decay  characteristics  of  Pu-238  and  Pu-239  are  listed  in  Table  5  and  provide  the  total 
energy  of  their  emissions.  The  predominant  mode  of  decay  for  Pu-238  is  with  a  5.49  MeV 
(average  energy)  a-particle  emission,  accompanied  by  low-energy  photons  and  electrons  (12% 
and  38%  per  disintegration,  respectively)  to  uranium-234.  Pu-239  decays  with  a  5.15  MeV 
(average  energy)  a-particle  emission,  accompanied  by  low-energy  photons  and  electrons  (4% 
and  32%  per  disintegration,  respectively)  to  uranium-235  (Shleien  1992).  Since  plutonium  will 
be  distributed  by  the  bloodstream  to  most  organs  throughout  the  body,  determining  how  much 
energy  is  absorbed  by  specific  organs  is  complex.  The  dose  delivered  to  each  given  organ  is  a 
result  of  energy  emitted  from  all  organs,  including  self-absorption  in  the  same  organ.  High- 
fidelity  phantoms  were  created  to  facilitate  calculating  these  doses  and  both  Monte  Carlo  and 
discrete  calculation  methods  were  implemented.  Oak  Ridge  National  Laboratory  (Christy  1987) 
generated  a  set  of  Specific  Absorbed  Fractions  (SAFs)  that  provides  an  energy-dependent  means 
of  calculating  the  dose  absorbed  by  an  organ,  as  a  function  of  the  energy  emitted  by  another 
organ.  For  example,  an  average  decay  of  Pu-238  emits  a  single  5.49  MeV  a-particle  (average 
energy).  In  addition,  Pu-238  decay  emits  a  photon  (average  energy  14  keV)  12%  of  the  time,  and 
38%  of  the  time  an  electron  (average  energy  22  keV).  Therefore,  for  every  Becquerel  (1 
disintegration  per  second)  of  activity  in  a  given  organ,  the  5.49  MeV  of  energy  from  the  alpha 
particle,  together  with  2  and  8  keV  (on  average)  of  photon  and  electron  energy,  are  emitted  per 
second. 


Table  5.  Decay  Properties  of  Pu-238  and  Pu-239  (NCRP  2008b) 


Half-life  (years) 

Average  Alpha 
Emission  per 
Decay  (MeV) 

Average  Photon 
Emission  per 
Decay  (MeV) 

Average 
Electron 
Emission  per 
Decay  (MeV) 

Pu-238 

87.7 

5.49 

0.002 

0.008 
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Pu-239 

24,110 

5.15 

0.001 

0.005 

By  interpolating  the  SAF  tables,  the  fraction  of  the  energy,  from  a  given  emission,  deposited 
in  another  (or  the  same)  organ  can  be  calculated.  The  contributions  from  all  organs  to  a  specific 
target  organ  are  summed.  For  Pu  isotopes,  the  RBM  and  lungs  are  considered  the  most  critical 
target  organs  for  predicting  acute  radiation  effects,  specifically  hematopoietic  syndrome  and 
pneumonitis  leading  to  lung  fibrosis.  The  equations  for  the  calculation  of  bone  marrow  and  lung 
doses  are  summarized  below  (note  these  are  generalized  equations  and  must  be  evaluated  for 
each  radiation  type  and  the  total  deposited  dose  summed): 
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where: 


DBM{t)  =  Dose  delivered  to  bone  marrow  at  time  t 
DL{t)  =  Dose  delivered  to  lung  at  time  t 
AlPu  =  Activity  of  Pu  in  organ  i 

Ej  =  Average  energy  emitted  per  decay  with  j  equal  to  the  radiation  types  1  (alpha),  2  (photon), 
or  3  (electron) 

At  =  Time  step  increment 

SAFjiPi  ->  BM )  =  Specific  Absorption  Fraction  from  organ  i  to  bone  marrow  with  j  equal  to  the 
radiation  types  1  (alpha),  2  (photon),  or  3  (electron) 

SAFj(Oi  ->  L)  =  Specific  Absorption  Fraction  from  organ  i  to  lungs  with  j  equal  to  the  radiation 
types  1  (alpha),  2  (photon),  or  3  (electron) 

Mbm  =  Mass  of  bone  marrow 

Ml  =  Mass  of  the  lungs 


The  SAF  values  for  the  low-energy  photons  and  electron  are  approximated  according  to  10 
keV  photons.  The  values  are  shown  in  Table  6  for  each  source  organ,  with  the  RBM  or  lungs  as 
the  target  organ.  The  electron  SAF  is  approximated  to  10  keV  photon,  since  it  is  the  closest 
relevant  SAF  available  for  use.  Some  assumptions  were  necessary  to  harmonize  the  SAF  values 
given  by  Christy  (1987),  with  the  Pu  tissue  content  provided  by  the  Leggett  (2005)  model.  For 
the  three  soft  tissue  compartments,  the  median  SAF  value  among  the  remaining  tissues  (skin, 
heart,  spleen,  etc.)  represents  soft  tissues  1,  2,  and  3. 


31 


Table  6.  Photon  SAF  Values  for  Pu-238/239  for  Red  Bone  Marrow  (RBM)  and  Lungs 


Source  Organ 

SAF  Values  -  RBM 

SAF  Values  -  Lungs 

Blood 

0.000447 

0.00193 

Lungs 

0.0000177 

0.913 

Liver 

0.0000291 

0.000132 

Kidneys 

0.0000154 

0 

Urinary  bladder  contents 

0.0033358 

0 

Cortical  bone 

0.0251 

0.00000944 

Trabecular  bone 

0.0639 

0.0000187 

Red  marrow 

0.0639 

0.0000187 

GI  contents 

0 

0 

Soft  tissue  1 

0.0000162 

0.000123 

Soft  tissue  2 

0.0000162 

0.000123 

Soft  tissue  3 

0.0000162 

0.000123 

Gonads 

0 

0 

Ovaries 

0 

0 

The  mean  free  path  of  an  a-particle  in  the  5  MeV  range  (as  is  the  case  with  both  Pu-238  and 
Pu-239)  is  measured  in  microns.  Thus,  it  is  assumed  that  the  alpha  particles  will  deposit  all  of 
their  energy  exclusively  within  the  organ  of  uptake.  The  RBM  is  a  complicated  case  because  of 
the  intricate  nature  of  the  structure,  in  and  around  the  blood  forming  region.  The  model  assumes 
the  a-particle  energy  is  deposited  locally  in  the  portion  of  bone  structure  where  the  Pu  is 
retained.  Therefore,  for  the  RBM  and  the  lungs,  it  is  assumed  that  all  of  the  alpha  particles 
deposit  all  of  their  energy  within  those  organs.  Likewise,  the  low-energy  electrons  have  a  very 
short  range  in  tissue  and  their  energy  is  assumed  to  be  absorbed  in  the  tissue  of  origin. 

4.5.2.  Effective  Whole  Body  Equivalent 

While  it  is  imperative  to  understand  doses  to  specific  organs  for  acute  effects,  most 
predictors  of  long-term  effects,  such  as  cancer  induction  from  internally  deposited  radioactive 
material,  are  based  on  the  concept  of  a  50-year  committed  whole-body  dose.  In  other  words,  the 
model  considers  the  radiation  dose  that  the  body  will  accumulate  over  a  50-year  period, 
following  the  uptake.  This  dose  value  is  based  on  the  estimated  dose  delivered  to  specific  organs 
or  tissues.  The  resultant  dose  is  then  multiplied  according  to  a  specified  tissue  weighting  factor 
(ICRP  1991).  While  our  model  approximates  the  calculation  of  dose  to  a  specific  organ  as  a 
function  of  activity  in  other  organs,  performing  this  calculation  from  all  organs,  to  all  organs,  and 
appropriately  weighting  the  results  is  not  practical,  SAF  values  are  not  available  for  all  tissue 
contributions  in  the  biokinetic  model.  A  large  number  of  assumptions  would  be  required  for  the 
calculation,  and  the  resulting  uncertainty  in  this  calculation  would  be  significant.  However,  we 
approximate  the  dose  using  tabulated  estimates  of  50-year  committed  doses  from  inhalation 
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uptakes  of  insoluble  forms  of  Pu-238  and  Pu-239.  By  using  some  basic  assumptions,  we 
developed  a  method  of  approximating  the  incremental  accumulation  of  whole-body  equivalent 
dose,  thereby  estimating  the  whole-body  effective  dose  over  time  and  the  50-year  committed 
dose. 


The  NCRP  estimates  a  50-year  effective  dose  equivalent  of  1.1E-05  Sv/Bq  (NCRP  2008b) 
for  inhalation  of  5  pm  insoluble  Pu-238  particulates  and  8.4E-06  Sv/Bq  for  inhalation  of  5  pm 
insoluble  Pu-239  particulates  (NCRP  2008b).  We  can  show  the  accumulation  of  this  dose  from 
the  following  equation: 


0(50) 


We  can  approximate  this  equation  with: 
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where: 


D(t)  =  Dose  rate  at  time  t 

If  we  assume  the  dose  rate  at  time  t  is  proportional  to  the  amount  of  Pu-238/239  in  the  body 
at  time  t,  then  we  can  rewrite  the  previous  equation  to: 
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where: 


APu(ti )  =  Plutonium  activity  retained  at  time  tt 
C  =  Constant  to  convert  APu  to  dose  rate 


Using  the  initial  retention  from  the  inhalation,  we  use  this  retention  equation  and  assume  the 
plutonium  delivers  a  “dose  equivalent  rate”  in  proportion  to  the  amount  of  Pu  retained  in  the 
body  at  time  t.  We  further  assume  the  Pu  is  in  equilibrium  throughout  the  body  at  all  times.  In 
other  words,  the  proportion  of  Pu  in  any  given  organ  is  constant  at  all  times.  We  then  calculate  a 
dose-rate  conversion  factor.  The  equilibrium  assumption  is  reasonable  for  the  majority  of  the  50- 
year  period  with  the  exception  of  the  period  immediately  following  the  uptake.  We  calculate  the 
proportionality  constant,  C,  by  numerically  integrating  the  results  from  our  model  of  Pu 
distribution  with  no  treatment.  Based  on  the  NCRP  equivalent  dose  estimates  for  the  50-year 
committed  dose  equivalent,  the  C  for  Pu-238  is  2.3E-09  Sievert/Becquerel/day  (Sv/Bq/day)  and 
Pu-239  is  1.75E-09  Sv/Bq/day. 

Using  this  conversion  factor,  we  can  estimate  the  accumulation  of  dose  equivalent  as  a 
function  of  the  Pu  retained  in  the  body  at  any  given  time,  regardless  of  the  elimination  rate.  This 
allows  us  to  compare  accumulated  doses  at  any  given  time  with  and  without  DTPA  treatment. 
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4.5.3.  Calculation  of  Acute  Health  Effects  Severity 

Acute  injury  may  be  possible  from  inhalation  exposure  to  plutonium  particulates.  Early  Pu 
workers  at  the  Mayak  Production  Association  in  the  former  Soviet  Union  were  exposed  to  high 
levels  of  Pu  through  inhalation,  and  radiation  pneumonitis/pulmonary  fibrosis  lead  to  many  of 
the  workers’  deaths  (Claycamp  2000).  Pneumonitis  and  pulmonary  fibrosis  are  considered  acute, 
deterministic  effects;  however,  the  time  scale  in  which  these  illnesses  manifest  are  longer  than 
many  other  acute  effects.  Furthermore,  the  onset  of  these  effects  is  dose  and  dose-rate  dependent. 
However,  the  precise  response  of  lung  tissue  to  Pu  deposited  in  the  lung  with  regard  to  dose  and 
dose-rate  is  not  well  characterized.  Therefore,  a  precise  model  for  acute  pulmonary  effects  is  not 
possible  due  to  limitations  in  data.  However,  a  severity  level  for  pulmonary  effects  can  be 
roughly  estimated.  Development  of  a  more  precise  injury  model  for  pulmonary  effects  is 
warranted. 

There  is  little  guidance  on  the  development  of  severity  levels  for  pulmonary  effects  as 
mentioned  previously.  Allied  Medical  Publication  8  (AMedP-8(C))  does  not  address  casualties 
from  internalized  radionuclides.  Cutaneous  and  external  whole-body-injury  severity  ranges  are 
not  applicable  to  the  effects  from  lung  doses  received  from  internalized  radionuclides;  however, 
whole-body  ranges  will  be  used  as  guidance  in  estimating  tissue-specific  dose  ranges  for  the 
severity  levels.  Mettler,  2008,  reviews  of  the  impact  of  radiation  on  all  organ  systems  provides 
limited  data  for  making  quantitative  assessments  on  predicting  incidence  of  pulmonary  effects. 
However,  such  data  is  needed  to  map  lung  dose  estimates  to  plausible  injury  severity  levels. 

Data  from  radiation  oncology  (Marks  2010)  and  from  occupational  exposures  (Claycamp  2000) 
were  evaluated  for  estimating  dose  ranges  for  pulmonary  tissue  effect  and  assigning  severity 
levels  1-3. 

Based  on  the  literature  reviewed,  the  symptomatic  incidence  in  50%  of  the  population  for 
pneumonitis  leading  to  fibrosis  is  assumed  to  be  approximately  8  Gy  prompt  exposure.  However, 
the  data  from  clinical  radiology  is  highly  variable,  greatly  dependent  on  dose  rate,  and  possibly 
has  confounding  factors  due  to  disease  and  treatment  (Marks  2010).  The  acute  radiation  dose 
which  causes  lethality  in  50%  of  the  population  (LD50)  for  the  RBM  is  2.9  Gy  (Anno  2003); 
accumulating  the  dose  over  a  longer  period  of  time  will  increase  this  value.  For  very  high  dose 
scenarios,  the  model  can  be  used  to  look  at  how  aggressive  and/or  effective  treatment  can  be. 

Based  on  the  literature  reviewed,  the  following  cumulative  dose  severity  levels  proposed  for 
pulmonary  effects  are: 

Severity  level  0:  0  -  <1  Gy 

Although  some  exposure  to  an  agent  or  effect  may  have  occurred,  no  observable  injury  (as 
would  be  indicated  by  manifested  symptoms)  has  developed.”  (AMedP-8(C)) 

Severity  level  1:  l-<5Gy 

Slightly  reduced  respiratory  capacity  in  a  portion  of  the  population 

Severity  level  2:  5  -  <8  Gy 

Transient  pneumonitis  resulting  in  some  respiratory  decline 

Severity  level  3:  >8  Gy 

Significant  pneumonitis  leading  to  fibrosis  and  eventual  pulmonary  failure 
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Note  that  the  proposed  ranges  have  considerable  uncertainty  associated  with  them,  since  data 
is  limited,  and  observed  effects  are  dependent  on  a  number  of  factors,  such  as  uniformity  of  dose, 
dose  rate,  etc.  A  separate  effort  specifically  aimed  at  establishing  reliable  ranges  for  pulmonary 
effects  from  inhaled  radionuclides  is  warranted.  Additional  uncertainty  can  arise  from 
extrapolating  the  anticipate  range  of  effects  from  prompt  doses  to  the  effects  anticipated  from 
radiation  delivered  over  a  long  period  of  time.  It  is  not  clear  whether  chronic  effects  of  fibrosis 
and  organ  impairment  would  occur  at  lower  or  higher  pulmonary  doses.  It  has  been  noted  that 
the  Mayak  worker  cumulative  pulmonary  doses  ranged  from  5-37  Gy  and  pneumosclerosis  was 
observed  in  several  participants  of  the  cohort  (Claycamp  2000). 

As  mentioned  above,  the  LD50  for  the  bone  marrow  is  about  2.9  Gy  (Anno  2003).  In  the  case 
of  inhaled  insoluble  Pu  radionuclides,  only  a  relatively  small  amount  of  Pu  is  absorbed  into  the 
systemic  circulation.  This  absorbed  Pu  would  become  fixed  in  the  bone,  but  the  total  bone 
marrow  dose  over  50  years  (from  a  1  pCi  intake)  would  be  less  than  1  mGy.  The  exposure  would 
need  to  be  significantly  larger  or  prolonged,  in  order  to  result  in  acute  effects  in  the 
hematopoietic  system.  Based  on  the  data  provided  by  this  modeling  effort,  acute  effects  to  the 
red  marrow  are  not  likely;  however,  increased  risk  to  long-term  health  effects,  such  as  cancer  is 
still  a  valid  concern.  Therefore,  the  dose  to  the  red  marrow  is  required  to  estimate  carcinogenic 
risk. 


4.5.4.  Calculation  of  Long-Term  Carcinogenic  Risk 


A  primary  concern  from  inhaled  and  internalized  radionuclides  is  the  increased  risk  to  long 
term  health  consequences  such  as  cancer  induction.  The  Biological  Effects  of  Ionizing  Radiation 
VII  (BEIR  VII)  carcinogenic  models  for  ionizing  radiation  are  the  current  standards  used  in 
several  countries,  including  the  United  States.  The  risk  model  developed  by  the  BEIR  VII 
committee  is  what  is  used  by  all  U.S.  government  agencies.  The  methodology  takes  into  account 
the  current  population  baseline  risk  to  cancer  and  projects  for  a  given  exposure  (radiation  dose), 
the  increased  risk  of  cancer  fatality  over  a  lifetime  (assumed  50-year).  However,  the  BEIR  VII 
risk  models  were  developed  for  external  ionizing  radiation.  To  apply  this  to  doses  received  over 
time  by  internalized  radionuclides,  additional  methodology  must  be  applied.  In  this  case,  the 
EPA  established  methodology  for  the  calculation  of  risk  (USEPA  1999).  Federal  Guidance 
Report  (FGR)  No.  13  establishes  a  method  for  calculating  the  dose  from  internalized 
radionuclides,  which  uses  the  integrated  dose  for  each  decade  with  the  age-respective  risk 
coefficients  provided  by  BEIR  VII.  Using  the  guidance  as  established  in  EPA  FGR  No.  13  and 
the  risk  models  from  BEIR  VII,  iterative  calculations  can  be  made  using  the  doses  received  over 
time,  with  and  without  treatment.  While  detailed  risk  analysis  can  be  conducted,  the  initial 
implementation  will  assume  the  exposure  occurs  for  a  20-year-old  healthy  male  or  female  (risk 
coefficients  vary  between  genders)  and  the  projected  risk  refers  to  the  increased  relative  risk  over 
50-years. 


IRR  = 


*RCi 


where  IRR  is  the  increased  relative  risk  to  cancer  mortality  over  50  years.  A  is  the  integrated 
dose  from  each  decade  after  exposure  using  the  effective  whole-body  dose  calculation,  described 
in  section  4.5.2,  and  obtaining  the  incremental  dose  for  each  of  5  decades.  RC  is  the  risk 
coefficient  for  the  corresponding  age  range. 
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4.6  Efficacy 


The  efficacy  of  DTPA  treatment  can  be  evaluated  in  terms  of  overall  dose  reduction  and  the 
reduction  of  risk  to  acute  and  long-term  health  effects.  A  dose  reduction  factor  is  calculated  by 
dividing  the  absorbed  dose  from  a  given  intake  of  Pu-238  or  Pu-239,  with  a  specified  treatment 
regimen,  by  the  absorbed  dose  from  same  intake  without  treatment: 

Dose  reduction  fCLCtOT  ^treated/ Duntreated 

where  Duntreated  is  the  radiation  dose  without  treatment  and  Dtreated  is  the  radiation  dose  with 
DTPA  treatment.  Note  that  the  two  doses  must  be  compared  at  the  same  time  post  exposure. 
Since  the  reduction  of  carcinogenic  risk  is  the  predominate  goal,  the  50-year  whole-body- 
committed  effective  dose  is  suggested  for  comparison  (see  section  5.2.7  for  further  discussion). 

The  efficacy  is  then  calculated  by: 

Efficacy  —  1  —  D  treated /^untreated 

where  Duntreated  is  the  radiation  dose  without  treatment  and  Ddtpa  is  the  radiation  dose  with 
treatment. 

It  is  important  to  note  that  DTPA  has  only  partial  efficacy  in  the  treatment  of  inhaled  PuC>2, 
due  to  its  limited  and  slow  absorption  of  PuC>2  into  the  systemic  circulation  from  the  lungs 
(Henge-Napoli  2000).  However,  PuCh  remains  the  most  available  form  of  plutonium  for  use  in  a 
radiological  weapon  and  understanding  the  fate  of  Pu  within  the  body  is  necessary  to  estimate 
radiation  dose  received  by  an  exposed  individual,  potential  health  effects,  and  long-term  health 
risks. 
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Section  5. 


Preliminary  Model  Results 

The  results  obtained  from  the  preliminary  implementation  of  the  model  are  presented  in  this 
section.  Final  results  may  vary  upon  final  implementation  in  stand-alone  code  and  software,  due 
to  the  precise  numerical  solutions  used  in  the  final  code  for  the  integrations  required  within  the 
various  models. 

5.1  Preliminary  Implementation 

A  number  of  the  calculations  required  in  the  Pu  decorporation  model  involve  several  sets  of 
differential  equations  and  the  integration  of  values  over  time.  Many  of  the  differential  equations 
are  interconnected  and  certain  calculations  require  simultaneous  numerical  solutions  to  provide 
the  appropriate  outputs.  Therefore,  Berkeley  Madonna™  v.  8.3.18  software,  licensed  by  Robert 
Macey  and  George  Oster  of  the  University  of  California,  was  used  to  solve  the  differential 
equations  contained  within  the  overall  model.  The  integration  method  used  was  the  2nd  order 
Runge-Kutta  method  and  a  time-step  interval  of  0.0001  days  was  employed  for  most 
calculations.  For  50-year  calculations,  a  smaller  time  step  of  0.05  days  was  used  for  optimal 
model  performance. 

Final  software  will  have  the  model  coded  autonomously  from  Berkeley  Madonna™  and 
incorporate  numerical  solutions  for  the  differential  equations  independently. 

5.2  Model  Outputs 

The  plutonium  decorporation  model  can  be  used  to  calculate  a  number  of  time  dependent 
parameters.  The  initial  conditions  currently  assumed  are  that  Pu  intake  is  through  inhalation  with 
a  given  air  concentration  of  Pu,  wind  speed,  breathing  rate,  and  average  particle  size.  The 
inhalation  and  uptake  model  requires  these  inputs.  If  values  are  not  known,  default  parameters 
may  be  used  for  wind  speed,  breathing  rate,  and  average  particle  size.  The  outputs  from  this 
portion  of  the  model  describe  the  levels  of  Pu  in  tissues  or  excretion  in  units,  to  use  as  the  input 
dose,  and  the  rate  is  per  day  (ex.  milicurie/day:  mCi/day): 

•  Amount  of  plutonium  deposition  in  the  three  major  areas  of  the  respiratory  tract 

•  Subsequent  plutonium  in  blood,  lungs,  GI  tract,  liver,  kidneys,  six  separate  bone 
compartments,  RBM,  urinary  bladder,  three  soft  tissue  compartments,  gonads  (or  ovaries) 

•  Total  plutonium  retained  in  the  body 

•  Plutonium  excretion  in  urine  and  feces 

•  Total  plutonium  excreted 

Certain  tissue  values  are  used  as  input  parameters  for  the  radiation  dose  calculations.  For 
radiation  dose  calculations,  the  input  dose  must  be  in  radioactive  units  (e.g.  pCi),  providing  an 
output  in  absorbed  radiation  dose  (pGy  or  pSv).  If  mass  is  used  as  the  input,  the  specific  activity 
of  the  exposure  material  can  be  used  to  convert  the  input  into  activity.  The  lung  and  RBM 
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calculations  are  derived  from  the  summation  of  the  high  and  low  LET  dose  contributions  from 
relevant  target  organs,  using  SAFs  as  described  in  section  4.5.1  and  integrated  over  the  time 
period  specified.  The  whole-body  effective  dose  is  calculated  based  on  the  decay-corrected  total 
Pu-238  or  Pu-239  activity  and  the  constant  for  the  whole-body  effective  dose  over  the  time 
period  specified,  as  described  in  section  4.5.2.  The  lung  and  RBM  acute  doses  are  used  to 
estimate  risk  of  acute  health  effects  (pulmonary  pneumonitis/fibrosis  and  hematopoietic 
syndrome)  and  the  whole  body  effective  dose  is  used  to  calculate  the  50-year  dose  for  evaluating 
long-term  carcinogenic  health  risk. 

Treatment  with  DTP  A  is  modeled  using  the  following  input  parameters:  dose  of  DTPA  (1 
gram,  fixed  parameter)  material,  which  results  in  a  vast  excess  of  DTPA  in  the  systemic 
circulation  as  compared  to  Pu,  treatment  start-time,  relative  to  the  initial  intake  of  Pu,  and  the 
duration  of  DTPA  treatment.  The  same  model  outputs  as  described  above  may  be  calculated  with 
adjustments  for  the  amount  of  Pu  chelated  by  DTPA  during  treatment. 

Untreated  and  treated  courses,  as  well  as  different  treatment  initiation  times  and  durations, 
can  be  evaluated  and  compared  to  determine  doses  and  dose  reductions  that  may  be  obtained  in 
different  scenarios.  The  data  generated  from  different  scenarios  are  detailed  below  to  illustrate 
the  reliability  and  utility  of  the  Pu  decorporation  model  developed  in  this  work. 

5.2.1.  Impact  of  Particle  Size  on  Plutonium  Biokinetics  and  Dose  Estimates 

As  discussed  in  Section  4.2.5,  the  particle  size  or  average  particle  size  distribution  of  the 
exposure  will  significantly  impact  the  observed  biokinetic  profile  after  an  intake.  Smaller  particle 
sizes  will  be  more  easily  inhaled  and  a  larger  fraction  of  that  exposure  will  reach  the  pulmonary 
region,  which  can  lead  to  systemic  uptake.  The  model  may  be  used  to  examine  the  impact  of 
particle  size  on  lung  and  systemic  body  burdens  and  clearance.  A  sample  evaluation  of  1  and  5 
pm  particle  sizes  with  two  different  breathing  rates  are  presented  in  Figure  1 1  and  Figure  12. 


Impact  of  Particle  Size  and  Breathing  Rate 


Figure  11.  Pu  Burden  in  the  Lung  as  a  Function  of  Particle  Size  and  Breathing  Rate 
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Impact  of  Particle  Size  and  Breathing  Rate 


Figure  12.  Pu  Retention  as  Function  of  Particle  Size  and  Breathing  Rate 


Since  particle  size  has  such  a  significant  impact  on  the  precise  deposition  and  resulting 
clearance  of  the  plutonium  intake,  the  particle  size  distribution  will  significantly  impact  radiation 
dose.  While  faster  clearance  with  smaller  particle  sizes  may  increase  the  systemic  body  burden 
of  Pu,  larger  particle  sizes  that  clear  more  slowly  will  result  in  a  higher  localized  exposure  to 
radiation  in  the  lung.  The  model  allows  one  to  examine  the  impact  of  particle  size  on  target  organ 
doses  as  well  as  whole-body  effective  doses.  The  impact  of  particle  size  (1  and  5  pm,  20  L/min 
breathing  rate)  on  radiation  dose  to  the  lung  from  a  1  pCi  intake  is  illustrated  in  Figure  13. 


Impact  of  Particle  Size  on  Lung  Doses 


1  micron  Pu-238 
1  micron  Pu-238 
5  micron  Pu-239 
5  micron  Pu-239 


Figure  13.  Radiation  Dose  to  the  Lung  from  Different  Particle  Sizes 
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5.2.2.  Plutonium  Tissue  Distribution 

The  model  allows  for  an  examination  of  the  biodistribution  of  Pu  throughout  the  body  over 
time.  Figure  14  illustrates  Pu  distribution  in  selected  tissues  over  30  days,  in  an  untreated 
scenario,  with  an  intake  of  1  pCi  Pu  (inhalation  exposure  to  PuCh,  5  pm  particles,  20  L/min 
breathing  rate).  Pu  has  an  affinity  for  the  liver  and  bone  compared  to  other  tissues.  This 
illustration  also  shows  how  Pu  has  an  initial  fast  clearance  from  the  respiratory  tract  for  a  portion 
of  the  intake,  which  then  slows  dramatically.  Pu  absorbed  from  the  lungs  redistributes  and  then 
accumulates  in  the  liver  and  bone. 


Biokinetic  Behavior  of  Plutonium  over  30  Days 


Time  (days) 


Figure  14.  Plutonium  Retention  in  Selected  Tissues  Over  30  Days 

The  total  Pu  in  Figure  14  is  calculated  by  summing  the  Pu  present  in  all  of  the  tissues 
represented  in  the  model  for  each  time  step.  Additional  organ  and  tissue  distributions  that  could 
not  be  easily  displayed  in  Figure  14  are  illustrated  in  Figure  15  and  their  distribution  is  extended 
out  to  180  days  to  show  the  cross-over  in  soft  tissue  distributions.  Figure  15  illustrates  that  the 
blood  has  a  minimal  amount  of  Pu,  having  been  redistributed  to  tissues,  which  retain  the  Pu  to 
varying  degrees.  The  different  soft  tissue  compartments  represent  tissues  that  have  rapid, 
intermediate,  and  slow  turnover. 
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Biokinetic  Behavior  of  Pu  in  Blood  and  Soft  Tissue 
Compartments 


Blood 
Kindeys 
Soft  tissue  1 
Soft  tissue  2 
Soft  Tissue  3 


Figure  15.  Pu  Distribution  for  Selected  Tissues  Over  180  Days 


5.2.3.  Comparison  of  Predicted  Plutonium  Retention  to  Experimental  and  Human  Data 

Data  collected  on  humans  indicates  large  variability  in  Pu  retention  and  biokinetics.  The 
precise  distribution  can  be  dramatically  influenced  by  the  chemical  form  of  Pu,  the  particle  size 
distribution  of  the  exposure,  and  breathing  rate.  Age,  gender,  and  body  weight  can  also  impact 
the  outcome. 

The  biokinetic  model  for  Pu  developed  by  Leggett  (2005)  and  implemented  in  our  model  has 
been  compared  to  experimental  animal  data  and  human  case  studies  (Leggett  2005,  NCRP 
2008b).  Figure  16  shows  Pu  retention  data  from  injection  studies  and  are  compared  to  the 
Leggett  2005  and  ICRP  1993  models  simulation  results. 


Time  after  injection  (days) 


Figure  16.  Model  Predictions  of  Pu  Retention  (Leggett  2005) 
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5.2.4.  DTPA  Decorporation  of  Plutonium 

DTPA  treatment  can  reduce  the  amount  of  Pu  in  the  body  by  binding  the  Pu  in  the 
circulation,  as  described  in  section  4.4.4.  DTPA  predominately  binds  Pu  in  the  systemic 
circulation.  Pu  is  then  re-equilibrated  according  to  the  variable  transfer  rates  of  the  different 
tissues.  The  impact  of  DTPA  on  Pu  whole-body  retention  is  illustrated  in  Figure  17.  This  figure 
shows  Pu  retention  over  60  days  after  1  pCi  lung  deposition  in  the  untreated  case  and  with 
DTPA  treatment  beginning  at  day  14  after  exposure,  with  standard  dosing  (1  g  daily)  for  a 
duration  of  30  days  (5  pm  particle  size,  20  L/min  breathing  rate). 

The  red  curve  shows  a  rapid  decline  of  Pu  after  treatment  is  initiated.  When  treatment  is 
discontinued  after  30  days,  the  slope  of  the  treated  curve  resumes  the  shape  of  the  untreated 
curve  and  is  parallel  to  it. 


Retention  of  Whole  Body  Pu  over  60  Days 


Figure  17.  Plutonium  Retention  With  and  Without  Treatment 


Large  individual  variability  is  observed  in  the  retention  and  excretion  of  Pu  with  DTPA 
treatment,  comparable  to  the  variability  observed  with  Pu  biokinetics  in  individuals.  Increased 
urinary  excretion  under  DTPA  treatment  ranges  from  5 -fold  to  as  much  as  150-fold,  depending 
on  the  exposure  conditions,  in  particular  the  chemical  form  of  Pu  involved  in  the  initial  exposure 
(Jech  1972,  Jolly  1972,  Schofield  1974,  Hall  1978,  Volf,  1978,  Wood  2000,  Menetrier  2005, 
Grappin  2007,  Bertelli  2010). 

5.2.5.  Plutonium  Excretion 

The  composite  model  estimates  cumulative  excretion  as  a  function  of  time.  The  excretion 
rate  of  Pu  changes  over  time  as  it  is  cleared  from  the  respiratory  tract  by  different  modes  with 
different  clearance  rates.  The  cumulative  excretion  over  60  days  from  the  scenario  presented  in 
the  previous  section  (5.2.4;  5  pm  particle  size  and  20  L/min  breathing  rate;  treatment  beginning 
day  14  and  lasting  for  30  days)  is  shown  in  Figure  18.  Total  excretion  is  dominated  by  fecal 
excretion  early  in  time  due  to  the  rapid  clearance  from  the  ET  and  TB  regions  via  the  GI  tract. 
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Excretion  of  Pu  in  Treated  and  Untreated 


Figure  18.  Plutonium  Excretion  With  and  Without  Treatment 

The  impact  on  urinary  excretion  during  treatment  for  this  scenario  is  shown  in  Figure  19.  The 
urinary  output  quickly  declines  after  treatment  ends.  The  output  after  treatment  is  stopped  is  less 
than  pre-treatment  values,  at  least  immediately  after  treatment  ceases,  due  to  the  time  required 
for  re-equilibration  from  slow  transfer  tissue  compartments  to  the  circulation. 


Urinary  Excretion  over  60  Days 


Figure  19.  Urinary  Excretion  of  Plutonium  With  and  Without  Treatment 
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Figure  20  compares  the  total  cumulative  excretion  to  fecal  excretion  in  the  treated  and 
untreated  case.  The  amount  of  Pu  cleared  through  the  GI  tract  is  significant,  as  illustrated  by  the 
magnitude  of  fecal  excretion.  Urinary  excretion  is  modest  without  DTPA  treatment  (initiated  on 
day  14  for  30  days),  but  dramatically  increases  upon  treatment,  as  evidenced  by  the  sharp  rise  in 
the  total  excretion  shown  in  the  graph  below.  The  rate  of  increase  excretion  quickly  returns  to 
that  of  pre-treatment  after  treatment  ceases.  A  slightly  lower  fecal  excretion  during  treatment  is 
due  to  the  decrease  in  systemic  plutonium  being  cleared  via  the  GI  tract  (and  being  excreted  in 
the  bile)  during  treatment. 


Pu  Total  and  Fecal  Excretion  over  60  Days 
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Figure  20.  Plutonium  Total  and  Fecal  Excretion  Over  60  Days 


5.2.6.  Plutonium  Acute  and  Whole  Body  Effective  Doses 

The  model  may  be  used  to  calculate  acute  doses  to  the  critical  target  organs  (RBM  and  lung) 
and  the  whole-body  effective  dose.  Since  both  Pu  isotopes  have  a-emissions,  high  radiation 
doses  delivered  directly  to  the  lung  tissue  is  possible.  The  RBM  is  one  of  the  more  sensitive 
organs  for  acute  effects  and  Pu  is  immobilized  in  bone  RBM  exposure  to  radionuclides  may 
cause  acute  affects,  such  as  the  hematopoietic  radiation  syndrome.  Our  model  can  determine  the 
dose  to  these  critical  organs  as  a  function  of  time  as  well  as  the  committed  effective  whole  body 
dose. 

The  acute  lung  and  red  marrow  doses  for  a  respiratory  deposition  of  lpCi  of  Pu-238  or  Pu- 
239  (5  micron  particles  with  a  breathing  rate  of  20  L/min)  in  untreated  cases  are  shown  in  Figure 
21.  The  acute  dose  refers  to  the  total  energy  deposited  in  a  particular  organ  and  potentially 
subsequent  acute  effects  (i.e.  cell  death,  which  results  in  organ  dysfunction  leading  to  morbidity 
and  mortality).  The  dose  to  the  lung  begins  to  decline  as  activity  is  cleared  from  the  region.  The 
dose  to  the  red  marrow  accumulates  as  Pu  is  redistributed  from  the  lung  and  incorporated  into  the 
bone  compartments.  Since  redistribution  is  slow,  the  red  marrow  dose  becomes  more  critical 
later  in  time. 


45 


Acute  Red  Marrow  Doses  after  l|iCi  Intake 
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Figure  21.  Acute  Doses  to  Critical  Organs  from  Pu-238  and  Pu-239  Over  90  Days 


The  concept  of  committed  effective  dose  was  developed  to  relate  long-term  health  effects 
such  as  carcinogenesis.  For  this  relationship,  tissues  must  be  weighted,  since  different  tissues 
have  different  dose  responses  to  carcinogenesis.  This  weighting  results  in  an  averaging  of  the 
dose  over  the  whole  body,  according  to  tissue-specific  sensitivities  and  provides  an  effective 
whole-body  dose  that  can  be  used  to  relate  to  overall  carcinogenic  risk. 

The  long-term  health  risk  to  cancer  is  most  often  evaluated  calculating  the  50-year 
committed  whole-body  effective  dose.  Table  7  shows  the  predicted  50-year  committed  whole- 
body  effective  dose  in  treated  and  untreated  scenarios  after  exposure  to  lpCi  of  Pu-238  and  Pu- 
239,  5  pm  particulates  (assuming  20  L/min  breathing  rate).  The  treated  scenarios  involve  a 
treatment  start  time  of  14  days  post-exposure  with  treatment  durations  of  60,  120,  and  180  days. 
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Table  7.  The  50- Year  Committed  Effective  Whole  Body  Radiation  Doses  from  Pu-238  and 

Pu-239  in  Untreated  and  Treated  Cases 


Dose 

Untreated 

DTPA  Day  14, 
60-day  Course 

DTPA  at  Day 
14, 120-day 
Course 

DTPA  at  Day 
14, 180-day 
Course 

Effective  Whole  Body,  Sv 
for  Pu-238 

0.29 

0.13 

0.11 

0.10 

Effective  Whole  Body,  Sv 
for  Pu-239 

0.26 

0.12 

0.10 

0.09 

The  data  in  Table  7  illustrates  reduction  in  the  50-year  committed  dose  by  approximately  55, 
62,  and  66%,  with  DTPA  treatment  beginning  day  14  and  continuing  for  60,  120,  or  180  days, 
respectively.  DTPA  treatments  afford  substantial  dose  reduction  from  Pu  radionculides,  allowing 
substantial  reduction  in  long-term  health  risks.  However,  given  limited  resources,  the 
calculations  may  help  optimize  treatment  regimens  for  maximum  reduction  in  Pu  body-burden 
and  dose.  In  this  scenario,  treating  an  individual  continuously  for  3  times  as  long  (180  days)  does 
not  result  in  three  times  the  reduction  in  dose.  In  fact,  the  need  for  re-equilibration  of  Pu  after  an 
initial  treatment  has  been  recognized  for  optimal  treatment  effectiveness. 

The  acute  dose  calculations  from  our  model  were  compared  to  the  current  approach  used  by 
the  ICRP  and  adopted  by  the  NCRP  (ICRP  1993,  NCRP  2008b).  Our  model  uses  the  same 
biokinetic  model  as  the  referenced  models;  however,  the  respiratory  clearance  model  differs  and 
is  specific  to  Pu  with  particle  size  dependence.  The  ICRP  and  NCRP  approaches  use  a  standard, 
generic  respiratory  clearance  model  developed  by  ICRP  (ICRP  1994).  Neither  of  these 
approaches  is  up-to-date  according  to  more  recent  recommendations.  For  Pu02,  the  ICRP  generic 
model  has  not  reliably  predicted  clearance  for  insoluble  Pu  radionuclides  (Bailey  2007,  Birchall 
2010).  The  model  implemented  in  this  work  uses  an  initial  clearance  of  Pu  greater  than  the  ICRP 
model  (as  recommended  by  Bailey  2007).  This  results  in  lower  lung  dose  and  greater  red  marrow 
dose  in  our  model,  as  compared  to  the  calculations  reported  by  the  NCRP  (2008b).  Even  small 
differences  in  respiratory  absorption  parameters  (shown  in  Table  2)  can  result  in  significant 
differences  in  estimated  radiation  doses.  The  range  of  values  referenced  in  the  literature  for 
respiratory  absorption  of  relatively  insoluble  plutonium  highlights  the  degree  of  uncertainty 
associated  with  current  plutonium  models. 
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Table  8.  Comparison  of  Doses  at  365  Days  (Gy/Bq  of  Pu-238  or  Pu-239  Intake)  by  Different 

Models 


Lung  -  Low 
LET 

Lung  -  High 
LET 

Red  Marrow  - 
Low  LET 

Red  Marrow  - 
High  LET 

NCRP  2008b 

Pu-238 

1.3x10-9 

1.5x10-6 

3.9x1 0-12 

1.6x10-9 

Model  Simulation 
Pu-238 

1.4x10-9 

8.1  x10-7 

1.3x10-10 

6.1x10-6 

NCRP  2008b 

Pu-239 

8.7x1 0-10 

1.3x10-6 

2.7x1 0-12 

1.5x10-9 

Model  Simulation 
Pu-239 

7.5x1 0-10 

7.6x1 0-7 

6.5x1 0"11 

5.7x10-8 

NCRP  model  simulations  were  based  on  inhalation  of  5  (am  particulate  for  workers  (not  stated,  but  assuming  ICRP 
breathing  rate  of  20  L/min  for  workers  was  used). 


5.2.7.  Evaluation  of  Efficacy 

The  efficacy  of  DTPA  treatment  can  be  evaluated  based  on  the  concept  of  dose  reduction. 
The  example  above,  in  Table  7,  is  a  good  illustration  of  how  the  comparison  of  doses  after 
DTPA  treatment  can  be  used  to  evaluate  efficacy.  This  type  of  analysis  can  allow  an  end-user  to 
determine  how  long  DTPA  treatment  administration  is  justified.  The  dose  reduction  from 
treatment  is  further  illustrated  in  Figure  22,  which  shows  the  dose  accumulated  over  90  days  in 
untreated  and  treated  cases  from  a  1  pCi  intake  of  Pu-238  or  Pu-239  (inhalation  exposure  to 
Pu02,  5  pm  particles,  20  L/min  breathing  rate).  The  treated  cases  assume  treatment  starting  day 
14  after  the  exposure  and  continuing  for  45  or  90  days.  A  significant  reduction  in  dose  is 
observed  with  treatment  for  45  days.  Dose  is  further  reduced  marginally  with  longer  treatment 
times  but  does  not  become  evident  until  longer  times  post-exposure.  Over  time,  as  the  dose 
continues  to  accumulate,  the  difference  becomes  more  apparent. 


Dose  Reduction  as  a  Function  of  Treatment  after  a 


Time  (days) 


Dose  Reduction  as  a  Function  of  Treatment  after  a 


-Pu-239 

untreated 


“  Pu-239  treated 
45  days 


“  Pu-239  treated 
90  days 


Figure  22.  Doses  from  Pu-238  and  Pu-239  Over  90  Days  in  Treated  and  Untreated  Cases 
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The  dose  saving  reduces  with  time,  which  illustrates  that  the  majority  of  the  dose  absorbed 
by  the  body  is  received  very  early  after  an  exposure.  The  time  of  treatment  with  DTPA  can  have 
a  significant  impact  on  the  overall  dose  resulting  from  Pu  exposure.  This  is  particularly 
important  for  actinides  such  as  Pu  and  Am  because  of  their  incorporation  into  bone.  After  Pu  is 
absorbed  into  the  systemic  circulation  from  the  lungs,  a  portion  can  quickly  be  retained  in  the 
bone  where  DTPA  has  minimal  ability  to  remove  it.  Having  DTPA  present  in  the  circulation  can 
help  bind  Pu  as  soon  as  it  is  absorbed  from  the  lung  and  reduces  the  amount  available  for 
incorporation  into  the  bone.  It  is  worth  noting  that  extended  and  intermittent  treatment  times  are 
often  necessary  due  to  the  slow  clearance  fractions  of  Pu  from  the  lungs.  Unfortunately,  DTPA 
injections  have  not  consistently  been  shown  to  effectively  increase  the  absorption  of  Pu  from  the 
lungs.  Therefore,  a  significant  portion  of  the  dose  observed  in  Figure  22  arises  from  the  Pu 
lodged  in  the  respiratory  tract. 

The  impact  of  the  different  delays  in  treatment  initiation  after  an  exposure  is  examined  in  the 
efficacy  of  decorporation.  In  the  example  shown  in  Table  9,  the  90-day  dose  from  1  pCi  of  Pu- 
238  is  reduced  more  dramatically  the  sooner  treatment  initiation  occurs.  Efficacy  of  treatment 
reduces  with  an  increase  in  treatment  initiation.  However,  because  DTPA  has  little  or  no 
influence  on  absorption  of  Pu  from  the  lungs,  no  significant  dose  reduction  is  afforded  to  that 
tissue  compartment.  Note  that  results  for  Pu-239  are  comparable  to  the  data  on  Pu-238. 


Table  9.  Doses  After  90  Days  from  lpCi  Pu-238  Exposure  in  Untreated  and  30-Day  DTPA 
Treated  Subjects  with  Different  Treatment  Initiation  Times 


Dose 

Untreated 

DTPA  at  day 

1, 30-day 
course 

DTPA  at  day 
14, 30-day 
course 

DTPA  at  day 
27,  30-day 
course 

DTPA  at  day 
40,  30-day 
course 

Lung,  mGy 

28.3 

28.3 

28.3 

28.3 

28.3 

Red  Marrow, 
mGy 

0.1 

0.02 

0.09 

0.12 

0.14 

Effective 
Whole  Body, 
mSv 

2.6 

1.3 

1.9 

2.3 

2.4 

Dose  reduction  is  further  examined  in  Table  10  where  the  50-year  effective  whole-body  dose 
versus  treatment  initiation  time  and  duration  of  treatment  are  presented. 
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Table  10. 50- Year  Committed  Whole  Body  Effective  Doses  (Sv) 


Duration  of 

Treatment  Initiation  Times  (days  post  exposure) 

DTPA,  days 

0 

30 

60 

90 

120 

150 

180 

0 

0.29 

- 

- 

- 

- 

- 

- 

30 

0.096 

0.22 

0.26 

0.27 

0.28 

0.28 

0.28 

60 

0.037 

0.20 

0.25 

0.26 

0.27 

0.27 

0.27 

90 

0.017 

0.18 

0.24 

0.25 

0.26 

0.26 

0.27 

120 

0.009 

0.18 

0.23 

0.25 

0.26 

0.26 

0.26 

150 

0.006 

0.17 

0.22 

0.24 

0.25 

0.25 

0.25 

180 

0.005 

0.17 

0.22 

0.24 

0.24 

0.25 

0.25 

The  data  illustrate  that  earlier  treatment  initiation  enables  greater  reduction  of  dose. 

Likewise,  longer  courses  of  DTPA  generally  provide  greater  dose  reduction.  However,  at  later 
treatment  initiation  times,  dose  reduction  is  not  as  significant  even  with  longer  treatment 
durations.  The  dose  reduction  factors  for  these  same  data  may  be  calculated  by  dividing  the 
treated  dose  by  the  untreated  dose.  For  dose  reduction  and  efficacy,  the  estimated  doses  of 
treated  and  untreated  cases  must  be  compared  at  the  same  time  post-exposure  for  a  scenario.  For 
all  the  scenarios  provided,  differences  in  the  50-year  committed  whole-body  effective  doses  are 
compared.  These  are  the  primary  endpoints  for  estimating  long-term  health  risk  (carcinogenesis), 
which  is  the  primary  concern  for  internalized  Pu. 

Dose  reduction  factor  =  D  D  oseufitrecited 

The  dose  reduction  factors  for  the  different  treatment  courses,  summarized  in  Table  10,  are 
listed  in  Table  11. 


Table  11.  Dose  Reduction  Factors  for  Different  Treatment  Regimens 


Duration  of 

Treatment  Initiation  Times  (days  post  exposure) 

DTPA,  days 

0 

30 

60 

90 

120 

150 

180 

30 

0.33 

0.76 

0.90 

0.93 

0.97 

0.97 

0.97 

60 

0.13 

0.69 

0.86 

0.90 

0.93 

0.93 

0.93 

90 

0.06 

0.62 

0.83 

0.86 

0.90 

0.90 

0.93 

120 

0.03 

0.62 

0.79 

0.86 

0.90 

0.90 

0.90 

150 

0.02 

0.59 

0.76 

0.83 

0.86 

0.86 

0.86 

180 

0.02 

0.59 

0.76 

0.83 

0.83 

0.86 

0.86 

The  values  in  Table  1 1  may  be  converted  to  estimates  of  efficacy  by  subtracting  the  dose 
reduction  factors  from  one. 

Efficacy  —  1  —  Dose  reduction  factor 
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The  efficacy  values  are  shown  in  Table  12.  The  higher  values  indicate  more  effective 
treatment  and  greater  efficacy.  Earlier  treatment  combined  with  longer  duration  affords  the 
greatest  dose  reduction  and  efficacy. 


Table  12.  Efficacy  of  Different  Treatment  Regimens 


Duration  of 

Treatment  Initiation  Times  (days  Dost  exposure) 

DTPA,  days 

0 

30 

60 

90 

120 

150 

180 

30 

0.67 

0.24 

0.10 

0.07 

0.03 

0.03 

0.03 

60 

0.87 

0.31 

0.14 

0.10 

0.07 

0.07 

0.07 

90 

0.94 

0.38 

0.17 

0.14 

0.10 

0.10 

0.07 

120 

0.97 

0.38 

0.21 

0.14 

0.10 

0.10 

0.10 

150 

0.98 

0.41 

0.24 

0.17 

0.14 

0.14 

0.14 

180 

0.98 

0.41 

0.24 

0.17 

0.17 

0.14 

0.14 

It  is  of  note  that  since  PuC>2  is  insoluble  and  is  lodged  in  the  respiratory  tract,  a  considerable 
portion  of  the  dose  from  inhaled  insoluble  radionuclides  of  Pu  cannot  be  alleviated  by  the  use  of 
DTPA  treatment.  The  use  of  DTPA  is  still  beneficial,  but  the  study  highlights  that  a  medical 
countermeasure  able  to  mobilize  and  remove  insoluble  Pu  radionuclides  from  the  respiratory 
tract  would  be  greatly  beneficial  in  the  reduction  of  dose  from  Pu  inhalation  exposures. 

5.2.8.  Acute  Injury  from  Plutonium 

Section  3.5.3  describes,  in  detail,  how  acute  injury  to  the  respiratory  region  is  possible  from 
inhalation  exposure  to  plutonium  particulates.  Estimates  made  from  our  composite  model  show 
intakes  of  Pu-23 8/239  in  the  mCi  range  would  be  required  to  develop  radiation 
pneumonitis/pulmonary  fibrosis  in  the  severity  ranges  2  or  3. 

5.2.9.  Long-term  Health  Effects 

At  lower  radiation  doses,  protection  from  long-term  health  effects,  such  as  carcinogenesis,  is 
the  primary  concern  of  internalized  radionuclides.  Each  branch  of  the  United  States  Department 
of  Defense  sets  guidance  and  recommended  dose  limits.  As  an  example,  the  U.S.  Navy  limits  the 
50-year  whole-body  committed  effective  dose  to  0.05  Sv  or  0.5  Sv  for  any  specific  organ  or 
tissue,  whichever  is  more  limiting  (NAYMED  P-5055).  The  model  can  help  to  resolve  the 
treatment  duration  required,  to  lower  the  committed  dose  within  these  limits.  The  recommended 
threshold  for  decorporation  treatment  is  any  exposure  resulting  in  a  50-year  whole  body 
committed  effective  dose  greater  than  200  mSv  (Rojas-Palma  2009).  Therefore,  the  model  can 
help  to  determine  whether  treatment  is  necessary. 

Using  this  approach,  the  increased  relative  risk  from  the  1  pCi  scenario,  presented  in  this 
section,  is  approximately  0.35%  for  Pu-238  and  0.28%  for  Pu-239.  In  the  case  of  prolonged 
treatment  (180  days),  the  risk  is  reduced  to  0.12%  and  0.09%  for  Pu-238  and  Pu-239, 
respectively. 
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Section  6. 


Limitations  of  the  Model 

The  composite  model  described  in  this  report  is  based  on  the  best,  currently  available  data; 
however,  as  with  all  models,  it  has  a  number  of  limitations.  The  purpose  of  this  section  is  to 
highlight  the  major  limitations  of  the  model  and  identify  possibilities  for  addressing  them  rather 
than  provide  a  comprehensive  review  of  uncertainties  associated  with  the  models. 

Large  variability  in  intake,  clearance,  and  biokinetics  after  Pu  exposures  has  been  observed 
among  individuals  as  well  as  in  animal  experiments.  Furthermore,  intake,  clearance,  and 
biokinetics  vary  significantly,  depending  on  the  specific  particle  size,  breathing  rate,  and 
chemical  form  of  the  exposure.  The  model  can  only  provide  a  tool  for  evaluating  specific 
scenarios  and  testing  different  treatment  regimens  with  fixed  assumptions.  For  an  actual 
exposure,  observed  bioassay  data  for  each  individual  must  be  used  to  determine  the  clearance, 
biokinetics,  and  resulting  dose  for  that  person. 

Application  of  the  inhalation  model  is  limited  to  wind  speeds  less  than  or  equal  to  9  m/s  and 
to  individuals  that  are  considered  normal  breathers.  The  biokinetic  model  is  currently  limited  by 
the  data  available  on  tissue  distribution  of  Pu.  It  assumes  three  generic  soft  tissue  compartments, 
which  clear  at  varying  rates.  If  actual  data  on  the  other  organs  were  available,  more  precise  dose 
estimates  would  be  possible.  More  detailed  data  on  specific  organ  burdens  may  become  available 
in  the  future,  as  autopsy  samples  from  the  United  States  Transuranium  and  Uranium  Registries 
(USTUR)  are  analyzed.  Based  on  animal  data,  it  is  known  that  a  fraction  of  the  less  soluble  Pu  in 
the  pulmonary  region  is  transferred  to  the  lymph  nodes  and  eventually  to  the  blood  stream  with  a 
long  half-life.  In  our  model,  the  lymph  nodes  are  not  treated  as  a  compartment  that  releases  Pu, 
due  to  a  lack  of  human  data  and  the  minor  contribution  of  this  transfer  mode. 

Finally,  our  model’s  removal  of  Pu  with  DTP  A  is  an  approximation,  based  on  the 
physiological  distribution  of  DTPA  and  Pu,  the  observed  data,  and  assumptions  concerning  the 
DTPA  binding.  The  approach  developed  in  this  work  affords  reliable  results  within  the 
limitations  and  uncertainties  of  any  of  the  existing  approaches  available  today.  Another  approach 
to  model  decorporation  with  DTPA  was  conducted  with  Pu  and  was  based  on  molecular 
chelation  kinetics  (Breustedt  2009,  2010).  In  this  approach,  the  molecular  ratios  of  Pu,  DTPA, 
and  their  binding  affinities  were  used  to  calculate  Pu-DTPA  formation  and  subsequent  excretion. 
The  investigators  reported  that  this  approach  did  not  provide  suitable  results  when  modeling 
repeat  DTPA  administration  over  a  short  time  period.  Also,  the  challenge  with  modeling 
molecular  chelation  kinetics  when  the  components  (Pu  and  DTPA)  are  present  on  vastly  different 
scales  was  discussed.  This  approach  may  still  afford  a  plausible  model  given  further 
development  and  experimental  data. 
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Section  7. 


Next  Steps 


Each  of  the  components  of  the  Pu  decorporation  model  will  be  consolidated  into  a  single  set 
of  codes  and  integrated  into  user-friendly  software  for  use  in  scenario  simulations.  These 
components  include  the  inhalation  exposure  model,  the  Pu  biokinetic  model,  the  DTPA  treatment 
model,  and  the  radiation  dose  models.  The  results  presented  here  provide  an  indication  of  the 
utility  of  the  composite  model  and  a  sample  of  the  type  of  data  that  can  be  generated. 

As  mentioned  previously,  large  individual  variability  was  observed  in  human  data  collected, 
to  date,  regarding  Pu  absorption,  retention,  and  response  to  treatment.  Although,  beyond  the 
scope  of  the  current  effort,  certain  aspects  of  individual  variability  could  be  addressed  in  the  Pu 
decorporation  model  presented  here.  For  example,  the  distribution  of  Pu  is  known  to  be 
dependent  on  age,  gender,  and  body  size/composition.  Distribution  is  partially  due  to  the 
differences  in  the  relative  size  of  the  tissues  and  partially  due  to  differences  in  metabolism. 
Physiological  data  exists  for  age  and  gender  specific  physiological  variances  (Richardson  2010 
and  ICRP  1995)  as  well  as  differences  in  long-term  effects  of  exposure  (BEIR  VII 2006). 
Extrapolations  can  be  made  mathematically  for  some  parameters,  to  account  for  body  size  and 
composition.  Therefore,  the  physiological  model  for  Pu  as  well  as  the  DTPA  treatment  model 
could  be  refined  to  account  for  physiological  differences  in  distribution  and  retention  of  Pu,  so 
that  additional  age  and  gender  specific  data  could  be  generated.  Likewise,  the  absorption  of 
radiation  is  impacted  by  the  differential  tissue  distribution  of  radioactivity.  Furthermore,  long¬ 
term  health  risks  vary  among  different  ages  and  between  genders.  Model  refinements  would  also 
enable  more  precise  health  risk  assessments  outside  of  the  “healthy,  adult  male”. 

As  mentioned  in  the  previous  section,  modeling  decorporation  of  Pu  on  the  level  of 
molecular  chelation  might  afford  a  more  precise  description  of  the  rate  of  removal  of  Pu  under 
DTPA  treatment.  However,  this  level  of  modeling  would  require  more  detailed  data  on  Pu  and 
DTPA  distribution  and  kinetics  in  tissues  than  is  currently  available.  The  effort  would  also 
involve  more  complex  modeling.  However,  with  additional  data  and  resources,  this  may  prove  a 
viable  approach  to  resolving  some  of  the  uncertainty  in  the  current  model. 
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Section  8. 


Conclusions 

The  DTPA  decorporation  model  presented  in  this  work  estimates  Pu  absorption,  distribution, 
retention,  excretion,  and  response  to  DTPA  treatment  in  adult  healthy  males  after  an  inhalation 
exposure.  The  model  is  easily  adapted  for  a  direct  ingestion  route  of  exposure,  which  may  be 
included  in  the  final  GUI.  The  model  further  estimates  the  acute  lung  and  RBM  dose  and  the 
whole-body  effective  dose  as  a  function  of  time. 

The  results  of  the  model  compare  favorably  with  both  human  data  and  alternative  models 
developed  by  the  ICRP  and  NCRP  (Figure  16  and  Table  8). 

Calculations  from  the  model  may  be  used  to  analyze  consequences  of  exposure  to  Pu-238 
and  Pu-239  (dioxide)  and  the  effect  of  treatment  based  on  initiation  and  duration  times.  The 
model  may  facilitate  interpretation  of  Pu  bioassay  data  and  aid  in  treatment  planning.  The  Pu 
decorporation  model  is  a  valuable  tool  for  assessing  the  effect  of  exposure  to  Pu-238/239,  in  the 
form  of  PuC>2  and  subsequent  DTPA  treatment. 

The  model  could  be  further  improved  by  incorporating  age  and  gender  specific  parameters 
and  data,  which  would  enable  more  precise  calculations  for  additional  segments  of  the 
population. 
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Section  10. 


Definitions,  Acronyms,  and  Abbreviations 


A1 

Alveolar  Interstitial 

AMedP-8 

Allied  Medical  Publication  8 

Am 

Americium 

ARA 

Applied  Research  Associates,  Inc. 

bb 

Bronchioles 

BB 

Bronchi 

BEIR 

Biological  Effects  of  Ionizing  Radiation 

Bq 

Becquerel 

CBD 

Chemical  and  Biological  Defense  Program 

CBI 

Information  Systems  Capability  Development  Division 

Ca 

Calcium 

CEDE 

Committed  effective  dose  equivalent 

Ci 

Curie 

Cm 

Curium 

Cs 

Cesium 

DFOA 

Desferrioxamine 

DoD 

Department  of  Defense 

DOE 

Department  of  Energy 

DTPA 

Diethylenetriaminepentaacetate  (or  diethylenetriaminepentaacetic  acid) 

DTRA 

Defense  Threat  Reduction  Agency 

EDTA 

Ethylenediamine  tetraaceticacid 

EEF 

Enhanced  excretion  factor 

EPA 

Environmental  Protection  Agency 

ET 

Extrathoracic 

EYS 

Environmental  Science  Division  (Argonne  National  Laboratory) 

Exp 

Exposure 

FDA 

Food  and  Drug  Administration 

FGR 

Federal  guidance  report 

g 

Gram 
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GI 

Gryphon 

GUI 

Gy 

ICRP 

IRF 

IUPAC 

HHS 

HP  AC 

JPEO-CBD 

JPM-IS 

JSTO 

KeV 

L 

LD50 

LET 

LN 

MCM 

MIRD 

mCi 

pCi 

min 

pm 

MeV 

MPPD 

NCRP 

NNSA 

NRC 

P 

Pu 

PUO2 

PBPK/PD 


Gastrointestinal  tract 
Gryphon  Scientific,  LLC 
Graphical  user  interface 
Gray 

International  Commission  on  Radiological  Protection 
Interstitial  removal  fraction 

International  Union  of  Pure  and  Applied  Chemistry 

Department  of  Health  and  Human  Services 

Hazard  Prediction  &  Assessment  Capability 

Joint  Program  Executive  Office  for  Chemical/Biological  Defense 

Joint  Program  Manager  Information  Systems 

Joint  Science  and  Technology  Office 

Kilo-electron  volt 

Liter 

Lethal  dose  in  50%  of  the  population 
Linear  energy  transfer 
Lymph  nodes 

Medical  Countermeasure  Model 

Medical  Internal  Radiation  Dose 

Millicurie 

Microcurie 

minute 

Micrometer,  or  micron 
Mega-electron  volts 
Multiple-Path  Particle  Dosimetry 

National  Council  on  Radiation  Protection  and  Measurements 

National  Nuclear  Security  Administration 

National  Research  Council 

Pulmonary 

Plutonium 

Plutonium  Dioxide 

Physiologically  Based  Pharmacokinetic/Pharmacodynamic 
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RBM 

Red  bone  marrow 

RD 

Research  and  Development  Directorate 

RDD 

Radiological  dispersal  device 

REMM 

Radiation  Emergency  Medical  Management  website  (HHS) 

RTG 

Radioisotope  thermoelectric  generators 

QF 

Quality  factor 

SAF 

Specific  absorbed  fraction 

SNS 

Strategic  National  Stockpile 

Sv 

Sievert 

TB 

Tracheobronchial 

TC 

Transfer  coefficient 

TERA 

Toxicology  Excellence  for  Risk  Assessment 

USTUR 

United  States  Transuranium  and  Uranium  Registries 

Zn 

Zinc 
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Appendix  A. 


Inhalation  Exposure  Model 


As  described  in  Figure  3,  the  inhalation  exposure  model  will  determine  the  amount  of 
plutonium  deposited  in  the  three  major  regions  of  the  respiratory  tract,  based  on  plutonium  air 
concentration,  duration  of  exposure,  wind  speed,  breathing  rate,  and  particle  size  distribution. 
Example  outputs  from  the  model  for  different  particle  sizes  and  breathing  rates  are  provided  in 
Table  13  and  Table  14.  In  Table  13,  the  fraction  of  the  exposure  that  actually  gets  deposited  in 
the  respiratory  tract  and  determines  the  intake  is  provided  for  two  different  particle  sizes  and 
breathing  rates.  In  these  cases,  intake  is  defined  as  the  amount  of  material  deposited  in  the 
respiratory  tract.  The  estimates  show  that  the  intake  is  significantly  impacted  by  the  particle  size 
of  the  exposure  and,  to  a  lesser  degree,  the  breathing  rate. 

Table  13.  Fraction  of  Exposure  Deposited  in  the  Three  Primary  Respiratory  Regions 


Region 

Particle  Size  and  Breathing  Rate 

1  pm  1  pm  5  pm  5  pm 

15  L/min  20  L/min  15  L/min  20  L/min 

Extra-thoracic 

0.006 

0.007 

0.208 

0.240 

T  racheobronchial 

0.061 

0.066 

0.266 

0.294 

Pulmonary 

0.093 

0.086 

0.260 

0.226 

Since  the  three  regions  clear  at  different  rates  and  by  different  modes  (see  section  4.2.5), 
particle  size  and  breathing  rate  further  impact  the  fate  of  plutonium  in  the  body  and  ultimately 
the  dose  received.  Larger  particle  sizes  and  faster  breathing  rates  yield  overall  greater  total  intake 
(see  Table  13). 

Calculating  the  intake  provides  the  fraction  of  the  intake  for  each  region,  (Table  14)  used  to 
estimate  the  amount  of  the  intake  in  each  of  the  three  respiratory  regions. 

Table  14.  Fraction  of  Intake  Deposited  in  the  Three  Primary  Respiratory  Regions 


Region 

Particle  Size  and  Breathing  Rate 

1  pm  1  pm  5  pm  5  pm 

15  L/min  20  L/min  15  L/min  20  L/min 

Extra-thoracic 

0.04 

0.04 

0.28 

0.31 

T  racheobronchial 

0.38 

0.42 

0.36 

0.39 

Pulmonary 

0.58 

0.54 

0.36 

0.30 
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Since  the  respiratory  tract  model  further  requires  additional  detail  about  the  distribution  of 
particles  in  the  tracheobronchial  region,  MPPD  was  used  to  determine  the  fractional  distribution 
of  particles  in  the  bronchi  and  bronchioles  in  the  TB  region. 

Table  15.  Fraction  of  the  Tracheobronchial  Deposition  in  the  Bronchi  and  Bronchioles 


Region 

Particle  Size  and  Breathing  Rate 

1  pm  1  pm  5  pm  5  pm 

15  L/min  20  L/min  15  L/min  20  L/min 

Bronchi  (TBbb) 

0.37 

0.44 

0.22 

0.27 

Bronchioles  (TBbb) 

0.63 

0.56 

0.78 

0.73 
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Appendix  B. 


Respiratory  Clearance  and  Uptake  Model 

The  respiratory  clearance  model  uses  the  initial  amount  of  plutonium  in  the  ET,  TB  (BB  and 
bb),  and  P  regions  to  calculate  the  amount  and  rate  of  transfer  from  the  respiratory  tract  to  either 
the  GI  tract  via  transport  processes,  such  as  ciliary  action,  or  directly  to  the  blood  stream  via 
absorption.  The  initial  amount  of  plutonium  in  each  of  the  respiratory  compartments  is  calculated 
by  adjusting  each  one’s  fraction  according  to  the  parameters  listed  in  Table  1  and  Table  13.  The 
values  for  the  TB  region  are  subdivided  into  to  bronchi  (BB)  or  bronchiole  (bb)  fractions 
according  to  the  partition  listed  in  Table  15. 

£T10  =  ETlfrac  *  ET 
ET20  =  ET2frac  *  ET 
ETseq0  =  ETseqfrac  *  ET 
BB10  —  BBlfrac  *  TBbb 
BB20  —  BB2frac  *  TBbb 
BBseq0  =  BBseqfrac  *  TBbb 
bbl0  -  bblfrac  *  TBbb 
bb20  =  bb2frac*TBbb 
bbseq0  =  bbseqfrac  *  TBbb 
AI10  —  AIlfrac  *  P 
AI20  =  AI2frac  *  P 
AI30  —  AI3frac  *  P 


The  movement  of  plutonium  particulates  within  the  respiratory  tract  is  depicted  in  Figure  6 
and  can  be  represented  by  the  following  mathematical  functions: 


—  (£Ti)  —  —ET1  *  TRET,^Env 


—  (ET2)  —  —ET2  *  TRET2-,gic  +  BBi  *  TRBbi-*et2  +  BBi  *  TRBb2^>et2  ~  ET2  *  fsr  *  sr  —  ET2  *  fss  *  ss 


d 

dt 


~JZ  (ETseq)  —  —ET  *  TRe Tseq^LNet 


dj.  (.BBi)  —  BB,  *  TRBB1^,ET2  +  bb,  *  TRbb,^BB,  +  bb2  *  TRbb2^,BB1  BB,  *  fsr  *  sr  BB,  *  fss  *  ss 
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The  compartments  are  defined  in  Table  1  of  the  main  text  and  respiratory  absorption  tract 
values  are  listed  in  Table  2  (Bailey  2007  values  are  implemented).  ET  is  the  extrathoracic,  TB  is 
the  tracheobronchial,  BB  is  the  bronchi,  bb  are  the  bronchioles,  AI  is  the  alveolar  interstitial,  and 
LN  are  lymph  nodes. 
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Appendix  C 


Plutonium  Biokinetic  Model 


The  plutonium  biokinetic  model  obtains  the  input  from  the  uptake  model  to  determine  the 
amount  of  plutonium  in  the  blood,  entering  the  systemic  circulation.  The  computation  is 
accomplished  by  solving  a  series  of  interdependent  differential  equations,  over  a  specified  period 
of  time: 


dt 


(PUbloodl) 


-Pu 


bloodl 


bloodl^liverl 


Pu 


bloodl 


bloodl^CS 


PU : 


bloodl 


*TQ 


bloodl^CV 


PU ; 


bloodl 


*  TCbloodl^>TS  —  Publoodl  *  PCbioodl^>TV  ~  Publoodl  *  TC bloodl^UBC  ~  P^ bloodl  *  TC bloodl^RT 
~  Publoodl  *  TCbioodl^OKT  ~  P ubloodl  *  P^bloodl^UGI  ~  Publoodl  *  PC bloodl^ gonads  ~  Publoodl 

*  PCbi00dl^ST2  ~  Publoodl  *  PCbi0odl^ST3  +  PUST1  *  PCsTl^bloodl  +  Publood2  *  PCbi0od2^bloodl 

+  ET2  *  fsr  *  sr  +  ET2  *  fss  *  +  BBi  *  fsr  *  Sr  +  BB1  *  fss  *  Ss  +  BB2  *  fsr  *  Sr  +  BB2  *  fss 

*  55  “h  bb 1  *  fsr  *  Sr  H-  bb 1  *  fss  *  Ss  *  fsr  *  Sr  ^^2  *  fss  *  *^s  AI^  *  fsr  *  Sr  AI^  *  fss 

*  55  +  AI2  *  fsr  *  sr  H-  ^4/2  *  ^55  *  ss  AI3  *  /sr  *  sr  ~h  AI3  *  /sS  *  5S  +  GIfrac  *  £"72  *  TCet2 -^gic 


(Publ00d2)  —  Publood2  *  PCbi0od2^UBC  Publood2  *  PCbi0od2^bloodl  Publood2  *  PCbi0od2^STl  +  PuOK 


P^liver3  *  TCliver3^>blood2  4"  P^gonads  *  P C gonads ->blood2  P^rM  *  TCt 


OK^blood2  +  PUST2  *  PCsT2^blood2  +  PUST3  *  PCsT3^blood2  +  Puliver2  *  P ^ liver 2^>blood2 

CM  *  1  u CM^blood2  +  PUTM 


*TC- 


TM^blood2 


fP^liverl)  PU-Uverl  *  P^liverl^SIC  P^liverl  *  P^liverl^liver2  4"  Publoodl  *  P^bloodl^liverl 

d_ 

—j^\PUuVer2J  ~  P^liver2  *  P C liver 2^blood2  P^liver2  *  PCiwer2->liver3  P^liverl  *  PC liver l^>liver2 

d_ 

~^\PUuver3 )  —  P^liver3  *  PCiiver3-*blood2  P^liver2  *  PCliver2->liver3 

d  , 

—  (PuST1)  —  PuST  1  *  7,C5ri^jbi00dl  +  PuMood 2  *  TCMoocL2^st1 
d 

—  (PuST2)  —  —PuST2  *  TCST2^hi00d 2  +  PuMood  1  *  PCbi00dl^ST2 
d 

—  (PuST3)  —  —PuST3  *  TCST3^bi00d 2  +  PuMood  1  *  PCbi00di^ST3 

d 

^(T^cm)  —  PUCS  *  PCcs^CM  +  *  PCcv^CM  ~  PUCM  *  PCcM^blood2 

d 

(Pucs)  —  ~PUCS  *  PCcs^CV  —  PUCS  *  PCcs-^CM  +  Publoodl  *  PCbloodl-*CS 

d 

—  (Pucv)  —  PuCs  *  PCcs^CV  —  PUCV  *  PCcv^CM  +  Publoodl  *  PCbloodl^CV 

d 

—  (PuTM)  —  PuTV  *  TCtv^tm  +  PuTS  *  TCts^tm  —  PuTM  *  TCTM^bi00d2 


—  (PuTS)  —  PuMoodl  *  TCbi00dl^TS  —  PuTS  *  TCts^tm  —  PuTS  *  TCts^tv 
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d 

dt 

d 

dt 

d 

dt 

d 

dt 


C PuTV )  =  PuTS  *  TCts^tv  —  PuTV  *  TCtv^tm  +  Publ00dl  *  TCbi00dl^TV 
(PiLrt)  =  PuMood i  *  TCMoocL1^rt  —  PuRT  *  TCrt^ubc 

( PuOKT )  =  Publoodl  *  TCbioodl->OKT  ~  P uOKT  *  P C0KT^blood2 
(■ P^gonads )  —  P^bloodl  *  PC  blooding  onads  ~  P^gonads  *  P C gonads ~^blood2 


df.  ( PUUBC )  —  PURT  *  TCrT^UBC  +  Publoodl  *  P^bloodl^UBC  +  Publood2  *  PCblood2^UBC  PUUBC  *  PC[JBC^urine 


d 

dt 


( PuUrine )  —  P^UBC  *  P^UBC^urine 


d 

—  (PuGIC)  —  PUnveri  *  PCu 


dt 

d 

dt 


liverl^SlC 


+  Pub 


bloodl  *  P^bloodl^UGI 


—  Plirnr  *  TC, 


GIC  *  7  ^GlC^feces 


+  ([PuET2T CET2^GIc)  *  (1  —  GIfrac ) 


(Pieces)  =  PUciC  *  ^ 


GlC^feces 


where: 


TC  is  the  transfer  coefficient;  ST1,  2,  and  3  are  soft  tissues  1,  2,  and  3;  CS  is  cortical  surface, 
CM  is  cortical  marrow,  CV  is  cortical  volume,  TS  is  trabecular  surface,  TV  is  trabecular  volume, 
TM  is  trabecular  bone  marrow,  and  RT  is  renal  tubules  and  OKT  is  other  kidney  tissue.  UBC  is 
urinary  bladder  contents,  UGI  is  upper  gastrointestinal  tract,  GIC  is  gastrointestinal  contents,  SIC 
is  small  intestines  contents,  and  GI/rac  is  the  fractional  of  absorption  of  Pu  from  the  GI  tract. 
Respiratory  parameters  are  defined  in  Appendix  A  and  B.  Values  for  transfer  coefficients  are 
listed  in  Table  3.  In  addition  to  the  Pu  specific  parameters  listed  in  Table  3,  the  generic  transfer 
coefficients  for  transfer  from  the  urinary  bladder  and  GI  tract  to  excreted  urine  and  feces 
(TCUBC^urine  and  TCCIC^feces)  are  7  and  0.5  per  day,  respectively  (ICRP  1993). 
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